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CHAPTER - 1 


INTRODUCTION AND BASIC 
CONCEPTS 

1.1 Introduction: 

In present study we consider Baysian analysis to software reliability 
models. The most important feature of using Baysian model is that prior 
information can be incorporated into the estimation procedure. Baysian 
methods are especially useful in reliability analysis because of the increase 
in reliability is usually achieved by improving a similar product developed 
previously. For software systems there is always some available - 
information about the software development and useful information can 
also be obtained from similar software products through for examples. 
Some measures of the software complexity. The fault removal history 
during the design phase etc. with an appropriate prior, Baysian inference 
is quite accurate and gives much better results than other methods such 
as method of maximum likelihood or least squares. Also Baysian method 
requires fewer test data to achieve a high accuracy of estimation. It shall 
be noted that care must be taken in choosing appropriate prior distributions 


(I) 


because an improper prior distribution may give very bad results which 
can even be worse than the maximum likelihood estimates. However, if 
no information seems to be available, non-informative priors can still be 
used. They usually give quite adequate results. However, existing models 
seem to be monotonic in a sense that simple and mathematically tractable 
posterior distributions are few in the existing literature. Lose of the existing 
Baysian models are also strangely related to the JM (Jelinski-Moranda) 
Model developed by Jelinski and Moranda ( 1 972), which is the simplest 
pure software reliability model and it seems quite popular, statisticians 
see the models discussed in Abdel-Ggaly et. al. ( 1 986) and Mazzuchi 
and Soyer (1988). It is argued that Baysian estimates of the parameters 
in the JM-Model usually have better behaviours than the maximum 
likelihood estimates and this is the reason why many such models are 
studied. Most existing Baysaan formulation of the JM-Model assume in 
a sense that the failure intensity per fault follows a Gamma distribution 
and the number of initial faults is Poison distributed. These assumptions 
are used mainly due to the analytical simplicity of the posterior 
distributions. There is a need to search for other suitable aprior 
distributions in order to describe other aprior information. A recent paper 
by Becker and Comarinopoulous ( 1 990) provides interesting ideas. Given 
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the values of x the time to the first failure is exponential with parameter 
X . The prior distribution ofhis assumed to belong to the following class 

n 

J {A) = c . It is a class of conjugate priors with some desirable 

properties, fhe original model is one which takes the possibility of a 
correct program into account and can be seem as a generalisation of the 
model by Thompson and Chelson (1980). Due to the computational 
complexity, it is important that posterior distributions are not too 
complicated. There are also some computational procedures developed 
see e.g. Mazzuchi and Soyer (1988). However as for the Langberg - 
Singhpurwallia's general Baysian model, the modelling procedures are 
essentially the same. The earliest paper presenting a Baysian Model is 
due to Littlewood and Verrall (1973). In the papers of Littlewood (1979a, 
1980 a) and (1980 b) where some other modifications may be found. 
Significant steps towards pure Baysian analysis of software reliability 
models have been taken by Singpurwallia and his colleagues Interesting 
papers are Meinhold and Singpurwallia ( 1 983), Craw and Singpurwallia 
(1984) Horigome et al (1984), Bahlow and Singpurwallia (1985) and 
Langberg and Singpurwallia (1 985). There are many papers dealing with 
Baysian treatment of the Jelinski - Moranda model. Important papers 
are Jewell ( 1 985), Littlewood and Sofer ( 1 987), Raftery ( 1 988) and Csenki 
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(1990). Other interesting papers Presenting or Discussing Baysian 
software reliability models are Thompson and Chelson (1980), Musa 
(1984), Kyparisis and Singpurwallia (1984), Ross (1985 a), Abdel - 
Ghaly et al (1986), Langberg and Singpurwallia (1 986), Jewell ( 1 986) Liu 
(1987), Mazzuchi and Soyer (1988), Catuneanu and Mihalache (1985), 
Horigome and Kaise (1990), Becker and CamarinopoTilous (1990) and 
Bunday and Al-Ayobui ( 1 990) 

1.2 A Brief Sketch of Statistical Techniques 
Used 

1.2.1 An Outline of Baysian Inference 

The foundation stone of this technique is 'Bayes' theorem and 
conditional probability. This theorem was first presented by Reverend 
Thomas Bayes (1963) and English minister. Laplace modified 'Bayes' 
basic theorem and the modified version is known as Bayes' theorem. 
Some of the standard texts giving more detailed discussions are those of 
Lindley (1965) and Box and Tiao (1973). For the basic theory and 
foundations one can refer to Sagave ( 1 962) and Jeffreys (1961 ). Bayes' 
analysis is an essentially self-contained paradigm for statistics. It is an 
excellent alternative to use of large-sample asymptotic statistical 
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procedures. Bayesian procedures are almost always equivalent to classical 
large sample procedure when the sample size is very large. 

In the parametric inference, the form of the population f(x,;i) is 
known while the parameter x known as labelling parameter, real valued 
or vector valued as the case may be, is unknown. However, we agree 
upon the parameter space i.e., the set of all possible values of the parameter 
which we denote by q . More so, this parameter x can also be represented 
by mean lifetime or failure rate or other such quantities associated with 
the underlying population from which the sample has been obtained. 

For example - the lifetime of a system follow exponential 

distribution with p.d.f. /( a-, A) = ye" 0,2 >0 (1-1) 

A 

X here stands for mean lifetime and the parameter space is q = 

{^U>o} 

In classical estimation theory the estimate of x depends only on 
sample values which we draw from f(x,A. ) and as such the information 
about X from other sources as well. If this additional information is in 
the from of a probability distribution of x ■, then it can be combined with 
data using Bayesian framework. There are cases in which x can be 


regarded as a r.v. with p.d.f. g(x ). For example in case of exponential 
model, the mean life x may be regard as varying from batch to batch 
over time and this variation may be represented by a probability 
distribution over q . As another example, suppose that woman has n 
sons, of whom x are haemophiliac (n-x) are normal. The probability of 
this event is 


p{x,X) 


( y{\ 


^Xj 


- xy-^ 


( 1 . 2 ) 


Where x is the probability that a particular sons will be 
haemophiliac. The problem is to make inferences about X. Given no 
additional information about a,, inferences would be based on (1.2) 
However, it may be possible to extract some information about x by 
examining the woman's family tree. For instance, suppose that the woman 
had normal parents but she had a brother who was hemophilic, then her 
mother must have been a carrier of the disease and therefore, she had a 
50% chance of inheriting the gene for hemophilic. If she did inherit the 
gene, then there is a 50% chance that a particular son inherit the disease 
( X =^ 2 ) and if she did not, all of her sons will be normal ( x =0). The prior 
probability distribution of x thus, becomes 


g(0) = p[2 = 0] = ^ 

g '( 2 ) = ~ 2^~2 


(1.3) 
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Thus the set up can be understood as follows: 

Suppose that n items are placed on a test. It is that their recorded 

life times form a random sample say x,, x^ , which follow a 

distribution with p.d.f f(x, To be specific we will assume x to be 
real valued. Consider x itself as a r.v. with p.d.f g( x )• Thus, the failure 
time p.d.f f(x,x,) can be regarded as a conditional p.d.f of x given x • 

i.e. f{x\X) where the marginal p.d.f of x is given by g(;^ ). Therefore 

the joint p.d.f of (x,, x^ x^, x ) is expressed as //(xjx, x„\A) 

= W' f{xi\X)g{X) = Z,(x,X 2 xJA) (1.4) 


The marginal p.d.f. of (x,, x, xj, is given by 


T(x,,a -2 .v„) = I //(x,,X2 x„\X)ci{A.) 

O. 


( 1 . 5 ) 


Thus the conditional p.d.f of x given data (x^, x, ..x^) is given by 


■>^2 



X„\X) 

T(x,,X2 X„) 


^(^^1.-^2 xJA)g(A) 

j L{xi,X 2 xJA)g(Z)dA 


( 1 . 6 ) 
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Thus, prior to obtaining the data (x,, xj, the variations in x 

were represented by g( x ) known as the prior distribution on x ■ However, 

after the data (x^, x^ x_^) has been observed, in the light of new 

information, the variations in x are represented by n(A.|x, ,X2 , the 

posterior distribution of x the uncertainty about the parameter a. prior 
to the experiment is by the prior p.d.f g(A,) and the same after the 

experiment is represented by the posterior p.d.f. / 

-(1.6) 

This process is a straight forward of Bayes theorem. After obtaining 
the posterior distributions of the parameters involved in the parent 
population, any statistical inference like estimation, testing of hypothesis 
about these parameters may be drawn with the help of this distribution. 

In Bayes approach, the estimator of x i.e., x is one which 

minimizes expected loss w.r.t. the posterior distribution i.e., it depends 

on the loss function chosen. If the loss function is taken as quadratic 

loss function defined as L (a , X ) = ix- X f then the Bayes estimator 

that accomplishes the task of estimating a. is the posterior mean, i.e., 

X* = E[X\x^,X2 JCJ = v„) dX 

n 

. " ' ■ ■ ■ ' ■ ' ' ■ ■ ■ ■ ' ■ 

A 


{H) 


Using the posterior distribution, a (1- a) 100% Bayes' confidence 
interval Pj, for x may be obtained from = 


cfA = I-( 


( 1 . 8 ) 


1.2.2 Concept ofPrior and Its Definition 

A detailed discussion to obtain the solution to the problem 
concerning the choice of a prior distribution of x is given in Raffia and 
Schlaifer (1961) but here we shall confine overselves just in defining 
them priors for the parameters in distributions differ in respect of their 
hidden properties. 


A family of priors is said to be conjugate, if it is 'closed under 
sampling [Lindly, 1972], If the same family, then this family is said to 
have 'closure under sampling' property [Weitherill, 1961]. Faiffa and 
Schlaifer, 1961] have considered a method of gathering prior densities 
on the parameter space. A family of such densities has been called by 
them a 'natural conjugate family'. For example in case of an exponential 
density, the inverted gamma priors from such a family. 




Jeffreys (1961) proposed a general rule for obtaining the prior 
distribution on the unknown parameter x . According to this rule, x 
which is assumed to be a r.v. follows such a distribution which is 
proportional to the square root of the Fisher information on x ■ 

Mathematically, org(x) = constant ^l(X) where, 

log i(;t|A.)|' 

1 I _ 

A difficulty arises when the prior information about the parameter 
is rather 'vague' or worse still there is no prior information whatever. 
This leads to the consideration of what are known as improper or quasi 
prior distribution. For a proper prior we have g( ;^ ) > o and j^g(A,)dk = 1 
while for an improper prior, g( A, ) > 0 but g(k)dk ^ 1 

Jeffreys prior may or may not be proper. Various rules have also 
been suggested for the selection of a prior but no neat solution appears 
to the problem till now. 


I{X) = E 


'aiogZ(;c|k)] 


dX 
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1.2.3 Non-Parametric or Distribution-Free Methods 

In life testing experiment it is not always practically possible to 
formulate a statistical model for failure time. Specially, at the early stage 
of designing a system, a designer knows little about the functional 
performance of the system. In such situations it becomes, important to 
develop distribution free techniques for the assessment of the reliability 
characteristics of a system. Kaplan and Meier (1958) suggested estimate 
for survival function based on an empirical study using censored data. 
Some other non-parametric methods have been discussed by Breslow 
and Crowley ( 1 974), Peterson ( 1 977), Johansen (1978) and Sharma and 
Krishna (1993). 

The most rewarding feature of non-parametric models is that we 
do not lose much in terms of answering the typical questions of concern 
in reliability studies, as compared to the more restrictive parametric models. 
Besides the obvious inherent advantage of the non-parametric methods, 
we also have the nice property of relative insensitivity of outliers in the 
data, superior power properties for a wide class of alternative distributions, 
and using the test statistics we can obtain estimators and distribution 
free confidence intervals for the population parameters of interest. 
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In fact, since F^, (x) is the sample mean of r.v's I(-c 3 o, x) (xi), ...I 
(-CO, x) (x^), we know by the central-limit theorem that F^(x) is 
asymptotically normally distributed with mean F(x) and variance 
— F(x)[l-F(x)]. Equation (1.10) and (1.1 1) show that for fixed x, F^(x) is 
an unbiased and mean-squared error consistent estimator of F(x) 
regardless of the form of F(.). 

1.3 Sequential Procedures 

The first fonnalised sequential procedure is the 'double sample' 
inspection plan by Dodge and Roming (1929) which was subsequently 
generalised to a multistage procedure by Bartky ( 1 943). Although there 
were other suggestion of multi-stage procedures (Hotelling, 1943), a 
systematic development of sequential methods started only in the last 
two years of World War II independently in USA and Great Britain. 
Abraham Wald developed a general theory of sequential analysis and in 
particular proposed the 'Sequential Probability Ratio Test' (SPRT) taking 
cue from a problem posed by the Statistical Research Group at Columbia 
University in 1943, Works of Wald (1947), Wald and Wolfwitz (1948), 
Arrow, Blackwell and Girshick (1949) and their co-workers laid the 
foundation of sequential decision procedures. 
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A procedure of making inferences about the distribution of one or 
more variable (r.v) is called a sequential procedure. The principal feature 
of such a procedure is a sampling scheme which lays down a rule under 
which one decides at each stage of the sampling whether to stop or to 
continue sampling, this decision being taken in the light of the observations 
already obtained. A general sequential procedure has two aspects first 
as a stopping rule, which tells us when to stop sampling and, second as 
an action rule which tells us what type of inference to make after sampling 
has been stopped. 

In order to judge the relative merits of two or more rival sequential 
tests, we make use of two criteria, one being the operating characteristic 
(OC) function and other the average sample number (ASN) function. 
The OC function often called power of sequential test, denote by L(0 ) 
is defined as the probability of accepting the null hypothesis when 0 is 
the true value of the parameter. It is closely related to the notion of the 
power function in the fixed sample size testing procedures. Now, since 
the number of observations say N, required by a sequential procedure to 
reach a decision is not predetermined but is a r.v. If we carry out the 
same sequential procedure repeatedly, we shall get different values of N. 


same sequential procedure repeatedly, we shall get different values of N. 
We may call a procedure preferable if it requires a small value of N on 
the average. This average value of N is called the average sample number 
of the sequential procedure. Thus, the OC function describes how well 
the procedure achieves objective of making correct decision, while ASN 
function represents the price one has to pay to reach a decision, in terms 
of the number of observations required by the test. 

1.3.1 The Sequential Probability Ratio Test (SPRT) 

Let X|, Xj .... be a sequence of identically and independently 
distributed (iid) random variables with common density (mass) function 
f(x,0). Let us consider the testing of a simple hypothesis 
against a simple alternative, H,:0 = 0, when the observations are taken 
sequentially. 

Let f^ denote the Joint p.d.fs (p.m.fs.) of x,, x, ...., x_^ under 

H,(i = 0,l) 

Define, ~ ''■’h) (1.12) 

Then, the SPRT for testing vs H, is a rule which states : 

(i) Accept Hq and terminate the experiment, if A,„(x) < B 



(ii) Reject and terminate the experiment, if A.„ (x) < B 

(iii) The experiment is continued by taking an additional observation 
if B<A,„(x)<A 

Here A and B (A > B) are the two positive constant which are 
determined so that test will have strength (a,p) where a,p are the 
probability of type I error and type II error respectively. If N is the 
stopping r.v., then 

a = (1-12) 

and p = P 0 ^(k„(x)< /4) (1-13) 

Here, A = ^ and B = — ^ 

a 1-a 


Efficiency of SPRT 


By the efficiency of the SPRT, say e, we mean the ratio 

= (114) 

' £ 0 ,.(.^) ^ ’ 


Here, N.(a,p) is the minimum value of E,,,( N) in the class (i=0,l ) 
clearly, 0 < e.< 1. An optimum test has'the efficiency one under H^ as 
well as under Hj. Wald (1947) has shown that the efficiency ofthe SPRT 
is very near to one under H^ as well as under H. 


1.3.2 Robustness 

Most of the statistical Inference problems are parametric in nature. 
We have assumed that the functional form of the distribution being sampled 
is known except for a finite number of parameters. It is be expected that 
any estimate or test of hypothesis concerning the unknown parameters 
constmcted on this assumption will petform better than the corresponding 
non-parametric procedure, provided that the underlying assumptions are 
satisfied. It is therefore of interest to know well the parametric optimal 
tests or estimates constructed for one population perform when the basic 
assumptions are modified. If we can construct tests or estimates that 
perform well for a variety of distributions, for example, there would be 
little point in using the corresponding nonparametric method unless the 
assumptions are seriously violated. 

In practice, one makes many assumptions in parametric inference, 
and any one or all of these may be violated. Thus one seldom has accurate 
knowledge about the tme underlying distribution. Similarly, the assumption 
of mutual independence or even indentical distribution may not hold. 
Any test or estimate that performs well under modifications of underlying 
assumptions is usually referred to as robust. 



The subject of robustness is receiving considerable attention of 
late. The first theoretic approach to the robust statistics was introduced 
by Huber (1964). His book on this topic Huber (1981) made this 
fundamental work to a wider audience. A different theoretic approach to 
robust statistics was originated by Hampel (1968, 1971, 1974). Robust 
statistics, in a loose, non-technical sense, is concerns with the fact that 
assumptions commonly made in statistics are at most approximations to 
reality. As a collection of related theories, robust statistics is the statistics 
of the approximate parametric methods. The moral is clear; One should 
check carefully to see that the underlying assumptions are satisfied before 
using parametric methods. 

1.3.3 Statistical Quality Control (Or Inspection) : A Way to 

Increase Reliability 

How can it happen that poor quality, unreliable products can be 
put out by a plant and then used by the consumer? In the first place, 
many of the properties of the unit that constitute its quality are either 
very difficult to check under factory conditions or are simply impossible 
to be checked because of the destructive nature of tests. In the second 
place the poor quality may be the result of the fact that the parameters 


characterizing the operation of the units are unstable, During the testing 
period, their values are completely acceptable to the user but then, after 
they stored or used for a while, they change in an erratic manner. 

Tests on durability and reliability are usually of destructive nature. 
However, inspection may be either destructive or “constructive. Since, 
testing the stability of parameters under different environmental conditions 
is of great significance in reliability theory, inspection is usually destructive 
or, in any case, it changes the quality of the units. These environmental 
conditions greatly affect the working of the engineering systems e.g., the 
debugging or bum-in period generally results in decreasing failure-rate 
and the wear-out period resulting due to aging-phenomena causes 
irreversible changes in the parameters or the unit leading to increasing 
failure rate. The failure - phenomena of an individual or a system is 
treated to be under normal operation when a constant or near constant 
failure-rate pattern is observed. After having an eye view of different 
failure-time models discussed earlier and in view of above, we can say 
that parameters in failure time distributions represent these environmental 
conditions. 
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Statistical Inspection (or S.Q.C.) problems are most closely related 
to the problem of ensuring reliability. Everything depends on the specific 
forms of the inspection operation. If a failure can occur during the 
inspection operation, S.Q.C. at the same time becomes a control over 
reliability. However, statistical inspection problems are more general. 
But properly restricted statistical inspection is a means of increasing the 
reliability of units, Obviously, during the course of inspections, testing 
techniques in statistical inference become the basis of decision making 
procedures on the disposal of allot of devices subject to constraints on 
risks (producer's and consumer's) and cost. 

1 .4 Posterior Analysis of the Mixture of Hazard-Rate Model 

Engineering Systems, while in operation, are always subject to 
environmental stresses and shocks which may or may not alter the failure 
pattern of the system. Suppose, p(0<p<l) is the unknown probability 
that the system usable to bear these stresses and its failure pattern remains 
unaffected, and also q(=l-p) is the probability that the system does not 
bear these stresses and the failure pattern of the system undergoes a 
change. In such situations, a failure distribution is generally used to 
describe mathematically the failure pattern of the system. To some extent. 
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the solution to the proposed problem is attempted through the mixture 
of distributions. To deal with such type of problem, Mendenhall and 
Hader (1958) have considered a failure population which may be divided 
into subpopulations, each subpopulation being exponentially distributed. 
The m.l.e's for the population parameter were obtained from samples 
censored at a predetermined test termination time. Kao (1959) siiggested 
the use of a mixture of two Weibull distribution as a lifetime model in 
which one of the components (Weibull) represents sudden (or 
catastrophic) failure and the other component (Weibull) represents 
wearout (or delayed) failure in electronic tubes. Proschan ( 1 963) explored 
the theoretical basis for a non-increasing failure rate and applied the 
results to real data from a mixture of exponentially distributed sub- 
populations. Kleyele and Dahiya (1975) studied a mixture of binomial 
and exponential population from censored data. The estimation of the 
parameters of a mixture of two distributions from a general parametric 
family of distributions assuming an ordering of the parameters was the 
subject matter of Shaked and Tran (1982). Papadopoulos and Padgett 
(1985) studied mixed failure population in which each subpopiilation 
was exponentially distributed MLE's as well as Bayes estimators for the 
population parameters and the reliability function were obtained from 


right censored samples. The Inverse Gaussina-Weibull mixture model 
was proposed as a life-time model by AL-Hussaini and Abd-El-Hakim 
(1989), in which each component represents a different type of failure 
such as the sudden or the wearout failure. It has been found that the use 
of the mixture model is essential in many situations to meet the relevant 
conclusions. It is mainly used to carry out analysis of weak models. The 
fertility analysis in survival studies is more adequate with the mixture 
model rather than the simple model. An interesting situation for using the 
mixture model is the analysis of there Reliability of a bulb knowing that 
the bulbs are being supplied by a fake company also with the original 
brand in the market. In this particular situation simple reliability can not 
help the analyst to reach a reliable conclusion. The additional advantage 
of the mixture model is that it helps in dealing with the more practical 
situations in the survival analysis. In most of the situations, while studying 
through a failure model we assume that the output is equally reliable for 
whole articles produced, while the guiding factor, i.e. machine or system 
may differ in context to produce the assumed articles. It is natural to say 
that all machines or systems in use may not produce the articles with 
equal reliableness and quality. Obviously, a new machine will produce a 
more reliable item than an older one. 


A lot of above discussion about the mixture model reveal that 
generally a failure distribution represents an attempt to describe 
mathematically, the life-time r. v, and failure data may reveal a number of 
physical causes responsible for the system failure. However, in this regard 
we face with two types of problems. Firstly, there are many physical 
causes that individually or collectively cause the failure of the system or 
device. Presently, it is not possible to differentiate between these physical 
causes and mathematically account for all of them, arid therefore, this in- 
isolation of the causes creates problem regarding the choice of a failure 
distribution. Secondly, even if a goodness of fit technique is applied to 
actual observations of time to failure, we face a problem a risint due to 
the non-symmetric nature of the life-time distributions whose behaviour 
is quite different at the tails where actual observations are appears inview 
of the limited sample size (Mann, Schaffer and Singpurwalla 1974). 
Obviously, the best one can do is to lookout for a concept which is 
useful for differentiating between different lifetime distributions. "Hazard- 
rate" which is the measure of instantaneous speed of failure, is one such 
concept in the literature on reliability. After analysing such physical 
considerations of the system, we can formulate a mixture of hazard-rate 
function which, in turn, provide the failure time distribution. 
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In view of the above, and due to continuous stresses and shocks 
on the system, let us suppose that the hazard-rate pattern of a with 
probability q. let the hazard-rate pattern of the system in these two 
situations be in either of the following two states. 

State 1 : Initially it experiences a constant hazard-rate (or chance 
failure) pattern and this pattern may (or may not) change with probability 
q(p=l-q). 

State 2 : If the stresses and shocks alter the hazard-rate pattern 
of the system with probability q. then it experiences a wear-out failure 
pattern. 

Using such a hazard-rate pattern, the characterisation of life-time 
distribution in the corresponding situations given. Various inferential 
properties of this life-time distribution along with the estimation of 
parameters and reliability characteristics is the subject matter of the study. 
Since, the estimates based on the operational data can be updated by 
incorporating past environmental experiences on the random variations 
in the life-time parameters [Martz and Waller 1 982]. Therefore, the Bayesian 


analysis of the parameters and other reliability characteristics is also 
given. Results have been highlighted with suitable example to show the 
application aspect of the developed procedures. 

1.5 Statistical Background 

Let, pc the probability of the event A, that the system is able to 
bear the stresses and shocks and its failure pattern remains unaltered. 

q = 1-p; the probability of the event complementary to A, i.e. x 
further, let, the mixture of the hazard-rate function be 
h(t) = pX + (1 - > 0,0 < p < 1 (1-15) 

In particularly, if, 

(i) p = 1, then the hazard-rate function h(t) = x and it represents 
the hazard-rate of an exponential distribution. 

(ii) p=0, then h(t) = xt and it represents the hazard-rate of the 
Rayleight distributioner Weibull distribution with its shape parameter as 2. 


1.5.1 Characterisation of the Life-Time Distribution 


Using, the well-known relationship ~ 


I 

I (x)dx 


and inview 


of (1.15), the p.d.f. of the life-time T can be obtained as 


f(t) = {pX + {l-p)lt) exp 


-J {/?X + (l-Tc)Xxlf£x 
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= [/?X + (l-p)?^r] + p)Xt^)-,k,t>Q (1.16) 

= 0, otherwise, 

the reliability function of the system is given by 

(1.17) 

and the mean time to system failure (MTSF) can be evaluated as 


R(t) = ^ = .exp[-{pXt + 1(1 - p)ll^ }]; t > 0 


MTSF = E(T) = J R(0dt = j exp[-{pXt + 1-(1 - p)Xt' }]dt 

0 0 

p^X 


= exp 


Jexpi 




2(1 -p) 

on making the transfonnation z = ^ 


l-pj 

m-p) 


dt 


(1.18) 


t + - 


P 


•P) 


in (1.18) and 


dz = A,(l - p)| 


t + - 


(1-p) 


dt , one gets 


MTSF = 


V2M1-P) 


exp 


2(1-P), 


r P-^ -l/2-U-r 


-10 


2(1 -p) 


dz 




expi 


p-X 


2(1 -p) 
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p -/. 

I -' V '”' -'-1 

To 1 


V ‘ 


0 


If "" ^ 

exp 

p-T.. 


'i-i 1 r 

(l2Ml-p)J 

_2(l-p)_ 


''-[2(l-p)J_ 


(1.19) 


where l.n(^)~ jz'”"‘e'''dz jg the well known incomplete gamma 
/ 

function and F 




= VTT 
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1.5.2 Classical Estimation of Parameters and Reliability 
Characteristics 

Let, t,, L , t be the random failure times of n items under test 

whose failure time distribution is as given in (1.16). Then, the 
likelihood function of the sample is given by 


= 


U{p + {l-p)ti\ 


i=\ 


exp. 




....( 121 ) 


1.5.3. 1 (MLE’s) Maximum Likelihood Estimates 

Using the likelihood function is (1.21) the estimates of various 
parameters involved in the p.d.f given in (1.21) can be obtained as 
follows 


Case I : When p is known 

Taking logarithm on both sides of (1.21), we have 


log/i'/,,/, S[p + (l-p)/,]-^Z ptf +-(l-p)^r ....(122) 

i~\ /=l[_ Z J 

Now, to obtain the m.l.e. of a, say £, we consider 

n 

....(1.23) 




51ogL(t|,t, ,t |?i.,p) n I • 1 • 

— ^ ^ = 0 which gives, ^ 

dX 


pti + ~^i\-p)h 


On using the invariance property of m.l.e's, the m.l.e. of R(t), h(t) 
and MTSF respectively given in ( 1 .22), ( 1 .26) and ( 1 .23) may be written 
as:- 

R(t), say (t), will be 

(1-24) 

(ii) The m.l.e. h(t), say h,(t), willbe h,(t) = i[p + (l-p)t] (1-25) 


(i) The m.l.e. for 


R,(t) = exp 


-X\ pt + -(l-p)t' 


(iii) 


The 


MTSF = 


71 


2X{\-p) 


expi 


m.l.e. 

2(1 -p) 


\-I 


for MTSF 

2{\-p) 


1 / 2 ' 


will be 
.,..,(1.26) 


Case II : When x is known 

To find the m.l.e. of p, say p, we consider 

aiogL(t,,t,....,t.,p.,p) 

dp 


(1-C) 


p+(l-p)l, 


■11 


h-2^F 


= 0 


orli 


T + 


0-X) 


-~lX(2-t,) = 0 

2 i=i 


(1.27) 


Estimate of p, i.e., p, can be obtained from (1.27) by using 
Newton-Raphson or other suitable iterative or inverse interpolation 


method. 


Again on using the invariance property of m.l.e's, one gets 


(i)Them.l.e. forR(t), say Rj(t),as R2(0-exp 


-X\ pt + 4(l-p)t 


..(1.28) 


(ii) The m.l.e. for h(t), say h, (t), as h, (t) = X[p + (l - p)t] 


(1.29) 


Case III : When x and p both are unknown 

Similarly, the simultaneous estimate of x and p are solutions of 
the equations 


\ = - 


or 


I 




(1-30) 


and 


Z 

/=i 


1 


A 


4zr,.(2-/,) = 0 


(1.31) 


(1.31) may be solved for p by Newton-Raphon or other suitable 
iterative method and this value is substituted in ( 1 .29) to obtain x ■ 


On using the invariance property of m.l.e.'s, 
(i) The m.l.e. for R(t), say R,( t) becomes 


(ii) The m.l.e. for h(t), say hjCt) , becomes 

h3(t) = £[p + (l-p)t] 


(1.33) 


1.5.4. Posterior Analysis 

Case I : When p is known 

Considering the conjugate prior distribution of x to be gamma 


withp.d.f. g 2 (^) = 7^^'^''e‘“';ot,p,X>0 (1-34) 

r(p) 

The likelihood function in (1.16) may be rewritten as 
Ltx,t2 = (1.35) 

Where, T,=n[p + (l-p)ti] (1.36) 

i= 


and T 


pti +^(1-P)t^ 


(1-37) 


Inview of the prior distribution of x in (1 .34) and the likelihood 

function in (1.35), the posterior distribution of x given t^, t, , t^^ may 

be derived as 


ni(^ihr^2 


T„p) 


7„|A,p)gi(^) 

j 

0 






)d/. 


0 


= ^ 2 )"!! x,''+P-ig->.(a.7; ) ; (a + 7, ),(« + P ), 7 > 0 

r(«+p) 


(1.38) 


Which is a gamma distribution with parameters (A. + ) and (n + P) . 

Therefore, the Bayes estimate of x , say x*, under the square- 
error loss function becomes, 


X' = ,t„p] = ,t„,p)dA, 

0 


(a + T)"^P 

e (a + T,)<a = 


n + P 
ot + T, 


(L.39) 


The corresponding posterior variance of x is given by 
V(X|t„t2 ,t„,p) = E(X'|t,,t, ,t,„p)-[E(X|t„t, ,t,„p)f (1.40) 


Now, E(A,"|(t|,t 2 t|,,p) — J ni(Nt|>t 3 ....»t 3 >p)d?'- 

0 


(O + r^^lj^n+P+le-Ma+T,)^^ - - + + P + 

r(n + P) I ' ~ (a + T,)' 


(1.41) 


Using (1.39) and (1.41) in (1.40), the posterior variance of x is 


V(X|t„t2 


4 nT) 


(n + P) 
(a + T)' 


(1.42) 


The Bayes estimate of R(t) is R'(t) = E[R(t)|t|.t, t„,p] 

CO 

= t„.p)dA 

0 

r(n+p) i 


f X 

1+- V 

a + I, 


yi+fi 


(1.43) 
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Where, % 


pt + (l-p)- 


(1.44) 


Also, = ni(A,|ri,i2 ,t„.P)d'R 

0 

1 


A 2T 
1 + -^ 




(X + T, 


(1.45) 


Therefore, the posterior variance of reliability function R(t) is given 
by V[R(t)|t„t 2 ,t„,p] = E[R'(t)|t„t, ,t„,p]-[E(R(t)it„t, ,t,„p)]- 

1 1 

\n+3 / „ \2(n+(5) 

(1-46) 


f 9T ^ 

1 + ^ 

V Ot + Tj y 


1 +- 




cx + T, 


2 J 


Similarly, the Bayes estimate of h(t), say, h’(t) , is 

h;(t) = E[t,,t2 ,t„,p] 

00 

= [p+(i-p)t]jxn,(xit„t,....,t„.p)dA 

0 

= [/? + (l-p)/].A, 

_ (n + P)[p + (l-p)t 
(a + T,) 


(1.47) 



and the posterior variance of h(t) is given by 

V[h(t)lt,,t,....,t„,p] = V[{p + (l-p)t}?A,42 t„,p] 

= [p + (l-p)t]^V(X|t,,t, ,t„,p) 

_ [p + (l-p)t]-(n + P) 

(a + T3)A 


.....(1.48) 
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Further once the posterior distribution of x has been obtained, 
one we get the Bayesian confidence limits for x and other reliability 
characteristics by making some suitable transformations. The interval 

(Pppj) is said to be a (1-a) 100% Bayes' confidence interval of xif 

00 

|ni(k|t|,t 2 ,t„,p)dX = l-a (1-49) 

Pi 

An equal tail (central) (1-a) 100% Bayes' confidence interval 
(PpP^), for X is given by 

Pi • “X 

Jn,(A.|t|,t2 ,t„,p)dk=-^= Jn,(y-itid2 ,t„,p)dk (1.50) 

0 ^ p. 

Class II : When x is Known 

Let, the prior distribution of p be Beta distribution with p.d.f 
giiP) = "x/ ' ^ * ~ <P<^ {1.51} 

1.6 Distribution-Free Inferences on System 
Reliability 

System reliability characteristics are analysed to achieve the 
preassigned reliability standards. As such, experiments are conducted 
to record failure time date which are used for analysing the systems in 
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terms of reliability function, increasing or decreasing hazard rate and 
mean time to system failure etc. Inferences on the reliability function are 
drawn in determining the worth of a system to accomplish and intended 
task. A vast literature dealing with this aspect is given in references like 
Kapur and Lamberson ( 1 977), Lawless ( 1 982), Martz and Waller ( 1 982), 
and Sinha (1986). In all such studies, an assumption about the form of 
the failure time distribution or the system or sub-system is made and it is 
observed that the reliability functions are usually complicated functions 
of the estimate of the failure time parameters. Also, at its early designing 
phase, it is not practically feasible to make assumptions about the failure 
time distribution of the system. More so, sometimes, it becomes 
necessary to form a system configuration in which the component may 
or may not be identical. In all such cases also it is difficult to state the 
lifetime distribution of the system. Thus, in many practical situations, it 
is iiicreasingly difficult to obtain failure time distribution of the system 
even if some admissible assumptions about the failure time distribution 
of the components are available. To overcome this difficulty, Sharma 
and Krishna (1995) obtained non-parametric distribution of the reliability 
function for some basic system configuration. The distribution so 
obtained may be used to get expected system reliability along with variance 


and higher moments. Further, since the life-testing experiments are usually 
costly and time consuming. Therefore, censored samples are used for 
the statistical treatment of data even at the cost of a loss in efficiency. In 
experiments using censored failure information, one would like to get 
sufficient information regarding the number of failures with a reasonable 
sample size or, time limit for drawing valid- inferences with a desirable 
accuracy. Here, one would like to restrict the sample size, while 
experimenting with high cost items, and limit the period of experimentation 
when the cost increases with time. The situation become more complex 
to draw parametric inferences about the reliability characteristics of the 
system when this time truncated failure information is used. It is therefore 
desirable to devise some procedure that are free of assumptions of the 
concerning failure time distribution. Some of the recent studies like Sharma 
and Bhutani (1992 b) and Sharma and Krishna (1994) have touched 
upon this aspect to some extent. Study in Sharma and Bhutani (1992 b) 
deals with the analysis of bounds for system reliability to perfomr an 
optimization between the sample size and censoring time. The 
developments in this study are quite useful for making an economic 
trade-off between sample size and censoring time. It is observed that 
estimates for the mean and variance of lifetimes distribution of the system 
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are enough to define such bounds. Still, the study has certain limitations 
as the bounds used for making such estimates of sample size and censoring 
time are restrictive to a situation when the mission time is greater than the 
mean lifetime. Sharma and Krishna ( 1 994) have used Kolmogrov-Smimov 
statistic for defining confidence limits for unreliability function. These 
confidence limits are found useful to define non-parametric estimates 
for the sample size and censoring time needed for the realisation of a 
preassigned failure information. 

Further, an engineering system may be better designed to meet the 
reliability objectives if the past environmental variations in the parameter 
involved in the distribution of a concerned variable are also incorporated 
with experimental data. In such a case, the non-parametric Bayesian 
analysis of the system reliability is of obvious importance. 

The development of some simple distribution- free techniques useflil 
for testing the significance of the system reliability for a given mission 
time or when a time truncated failure information is to be used. The test 
fortesting the significance difference between the reliability functions of 
the two systems with common mission time and for testing the significance 


of a specified difference between two reliability functions are also given. 
The Bayesian estimation of the system reliability with a suitable prior 
distribution has been discussed here. 

1.7 Bayes Test Plans for Testing A Series 
Configuration 

A good manufacturer for the acceptance purposes always keeps 
a constant watch on the quality of products he manufactures. Such 
alertness enables to specification of the purchaser. On the other hand a 
purchaser is also anxious to satisfy himself about the quality of products 
he accepts. In quality control texts the manufactures and the purchaser 
are referred to as producer and consumer respectively. In order to assess 
the quality of products, the acceptance inspection is usually based on 
sampling. All acceptance plans for the items which are destructive in 
nature must inevitably be done by sampling. In many other instances 
sampling inspection is used because the cost of 100% inspection is 
prohibitive. Moreover, when there are a great many similar items of 
product to be inspected. It has been found that if a scientifically designed 
sampling inspection plan provides protection to consumer against the 
acceptance of too much defective product and, on the other hand. 
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producer is protected against the rejection of too much good product. It 
is obvious that consumer and producer have completely opposite view 
points towards the selection of sampling plans, i.e., they have conflicing 
interests. Practically it has been observed that no sampling plan can give 
complete protection to producer and as well as consumer, i.e, a sampling 
inspection plan is likely to involve some amount of sampling error which 
with result in producing certain risks on the part of producer and consumer 
both. These risks will be needful in construction of an acceptance 
sampling plan. Further, variation in the quality of manufactured product 
in the repetitive process in industry is inevitable. This variation will intern 
affect the producer's and consumer's risks. Considering this fact. Brush 
(1986) presented the process of comparing Bayes acceptance sampling 
plans with classical sampling plans. The study highlighted the importance 
of using Bayes producer's risk (BPR) as a supplement to the analysis of 
the modified classical measure of producer's risk (CPR*). After 
normalizing the quality on the index scale as in Brush et. al ( 198 1 ) and 
Hoadly (1981), Brush (1986) defined BPR and compared it with the 
CPR*. The study in Brush (1986), however, does not define Bayes 
consumer's risk (BCR) and the modified classical measure of consumer's 
risk (CCR*). The difficulty in defining them concerns with the appropriate 
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definition of bad quality needed for defining BCR and CCR*. Later, 
Sharma and Bhutani (1992 c) investigated an acceptance sampling plan 
with given producer's and consumer's quality specifications. In addition, 
the concept of BCR is analysed as a supplement to CCR*. 

Further, various engineering systems are analysed in respect of 
their reliability characteristics to meet the established reliability standards. 
It has been felt that the reliability characteristics of the system can be 
better analysed if the past environmental variations in the lifetime 
distribution of the system are also in the analysis. Similar to product 
control problems in the statistical quality control (SOC), systems are 
tested for their reliability characteristics like. MTSF, reliability, availability 
and hazard rate etc. Likewise, acceptance sampling plans need to be 
developed for accepting or rejecting a lot of system on the basis of 
certain reliability specifications when the parameters of the lifetime 
distribution are also taken into consideration. 

Following the above concepts, the present chapter deals wit the 
development of a Bayesian life testing plan used for testing the lots of 
items or subsystems that are to be arranged in a series configuration. 
The and producer's risks. The concepts have been highlighted and 
analysed with example. 


( 39 ) 


CHAPTER - 2 


BAYESIAN ANALYSIS OF FIRST 
CONCEPITION TO LIVE BIRTH 

2.0 INTRODUCTION 

At the time of marriage a woman is susceptible to conception and 
the time elapsed before a conception is a random variable detennined by 
fecundability, which is define as the monthly chance of a conception. It 
is important here to note the time of first conception after the marriage 
because the analysis of waiting time of first conception signifies couple's 
fertility at early stages of married life. This variable is widely used to 
study fertility characteristic of a woman, since it is independent of effect 
of amenorrhea period and generally, a woman does not like to use 
contraceptives to postpone first birth. There is little chance of recall 
lapse in reporting the time of first birth from the date of marriage of first 
birth. 

Treating the first conception as a random phenomenon the 
probabilistic models can be developed. And this variable can be treated 
as discrete as well as continuous variable depending upon the situation 

and assumptions made for the study. 
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For the first time this variable was considered as a discrete variable 
and Gini ( 1 924) derived the geometric distribution for the waiting time of 
first conception. He defined the term 'fecundability' as the monthly chance 
of conception for woman living in the married, fecund and exposed 
state. 

Later the same variable was considered as a continuous variable 
and hence treating the same elapse from the marriage or from the begimiing 
of the reproductive process to first conception as continuous makes 
mathematical treatments more convenient and easy. Singh ( 1 964), Henry 
(1953) and Vincnet (1961) developed some models treating the waiting 
time of the first conception as continuous. The negative exponential 
distribution plays the role of geometric distribution for studying the waiting 
time of conceive after marriage. Thus, if T denotes the time of first 
conception then its density function, say, fft) is given by 
f(t) = 5e'"’;t > 0,5 > 0 Where 5 is instantaneous fecundability. 

A number of authors made modifications on the above simple 
distribution to study realistic situations. 

In the present model of waiting time of first conception the time 
elapsed is defined over the range (0, co). But in practical problems the 
upper limit may be considered as finite; as a woman can conceive upto 
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an age limit. So, there is a need of introducing a new continuous model 
with finite range. Keeping this in view, an attempt has been made to 
characterize an existing model derived by Mukeherjee Islam ( 1 983 ), defined 
over a finite range for the puipose of life testing analysis but it suits in 
realistic or real life situations. 


2.1 THE MODEL 


A new probability distribution has been considered in this section 
as a continuous model, introduced by Mukhejee-Islam (1983) for the 
purpose of studying waiting time. 


f(t;0,p) = 


t'"'; p,e>o; 


t>0 


The above model is monotonic decreasing and highly skewed to 
the right. The graph is J-shaped thereby showing the unimodel feature. 
The same model has been modified on the range of parameter '0' 


for the specific use in the study. 


f(t;0.p) = 


f_Ll 

Uv 


-p-i 


G>ri; 0<t<(); 0<p<l 


( 2 . 1 ) 


Where 'p' is instantaneous fecundability and '0' is considered as 
age limit beyond which a married woman can not conceive. Here w'e 
define the 'p' over the range (0, i ) as it has been defined in thd present 
model as instantaneous fecundability. 
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The distribution function of above model will be 


F(t) = P(T < t) 

t 

= jf(t)dt 
0 



lev 


p 


Jo 


or finally 


The survival function at time t, say S(t) - P(T>t) is given by 


I 



S(t) = p[T > t] 


= l-P[T<t] 


by 


= 1 - F(t) 

^l_riT (2.3) 
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Also, the conception rate function, say, w(f) at time t is then gi\ en 


\v(t) = 


f(t;6.p) 

P(T>t) 
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Which finally comes out to be, 


W(t) = 


Pt 


p-i 


0P_t'’ 


.(2.4) 


2.2 Characterization of The Model 

The important main features of the proposed model are found as 
follows; 

2.2.1 Mean Time of First Conception 


u 

E(T) = Jt.f(t)dt 



t| f’-'dt 
. 6 ”, 

V0P 


_P_ 

0 '’ 


. p+i 


p + f:., 


i.e. E(T) - 


P- 


.(2.5) 


2.2.2 VARIANCE OFTIME OF FIRST CONCEPTION 

In order to calculate Itratly we calculate lc{T’ ) as; 

B 

E(T-) = Jt-f(t;0,p) dl 


0 

e 


dt 


0P 
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.(2.6) 


which finally yields; 

E(T-)=^e- 

p+1 

Putting the values from (2.5) and (2.6) in the following expression; 
V(T) = E(T-)-[E(t)]' 


P -0^ 


p + 2 


P + 1 
• 0 - 


(P + l)(P + 2) 


.(2.7) 


2.2.3 MAXIMUM LIKELIHOOD ESTIMATE 

The likelihood function for the model (2. 1 ) is given by 

L(t;0) = p''0'''»’flt|’^' 

1 = 1 

LogL{t;0) = n logp- nplog0 + (p - 1)^ logt, (2.8) 

Differentiating the above equation partially with respect to 'p' and 

equating it to zero, 

5LogL(t;0) n , ,, Vi , 

ii— = — n lot>0 -f > loa t. = 0 

<30 p - A. - : 

The m.l.e of p is finally obtained as 

‘’“nloge-Vklgi; 

Again, differentiating the equation partially with respect lo '()' and 

equating it to zero to obtain the m.l.e of o 
(3LogL(t;0) np _ 


W 
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In the solution for MLE of 0 the traditional method is not applicable. 
The MLE is obtained through order statistic technique. Since the upper 
limit of the model is 0 , it is convincing to take t^^^^ i.e. maximum t. as the 
m.l.e for the parameter 0 

i.e. 9 = = Max(t,,t 2 ,t„) (2. 1 0) 


2.2.4. MEDIAN OF THE TIME OF FIRST CONCEPTION 


By the definition of Median, we have 
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Hence, 



2 ’’ 


( 2 . 11 ) 
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2.2.5 r " MOMENT ABOUT ORIGIN OF TIME 
OR FIRST CONCEPTION 



2.3 BAYESI AN ANALYSIS OF THE MODEL 

In the present section an attempt is made to fiild out posterior 
distribution of the model discussed in previous section. In that particular 
model, there are two parameters p and 0 , where p stands for fecundability 
and 0 stands for upper age limit. Let us suppose that q varies from 
woman to woman according to social, economical, biological, 
enx'ironmental etc. factors, where p the fecundability is constant. 

The three prior distributions to be considered for the purpose ol' 
the Bayesian analysis are Uniform, Inverted (iamma and the one 
proposed by Siddiqui et. al belongs to Unite range of family of 
distributions while inverted Gamma belongs to infinite range of family of 
distributions. In the past literature it was observed that the results presented 
were complicated and not get very handy to solve out numerically. The 
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finite range priors are successfully tried in this article in order to reduce 
the complexity of the results. 


are 


The probability density function of the model is discussed i.e. 
M.I. model is as follows 

f(t;e,p) = (p!e'’)t'’-‘ (2.13) 

and the suiwival function and conception rate function of the model 

S(t) = l-(tj0)' (2.14) 

and 

w(t) = -^ (2.15) 

respectively. 


2.3.1. UNIFORM PRIOR DENSITY 

Now let us suppose that 0 follows Uniterm distribution. 1 he 
probability density function ofUnifonu distributions; in the interval [(x.p] 

truncated at 'a' is gi\ en by 

(a -!)(«- (I)""' 1 


2 ( 0 ) 


, ; 0<a<B<a 


( 2 . 16 ) 


0 ^' 

Now, suppose t|, U , t„ denotes n observation then the likelihood 

conditional function of 0 is given by 


(48) 


L(eitpt2 t„) = nf(t,|0) 




p-l 


!) 


= p''e'"pntr 

= H.p'’0“ 


i=l 
.iA~np 


.( 2 . 17 ) 


where H=ntf' 

i=l 

The likelihood conditional function under the assumption of prior 
density (2. 1 6) provides the following posterior distribution 


G(0|t|d2 d,,)^ 


H.p'’0 


nn-np(ayJ)(aP>"“' 




J 0^' 

a ' 


-up- a 


je^"'’"’d0 


Which solves out to be in the follow ing form 


G(0|t|,t, tj = 






(I - a - up) 


.\s 


N. 


where 


( np + a - 1 ) 

T 

(V’’’ ^ 


.( 2 . 18 ) 


m 


The Bayes estiiTiate of the parameter 6 can be obtained as 


e- = feide 

a 


N, 

' 1 1 ' 

^np+a-2 p np+a-2 

(np + a-2) 


(np + a-l(ap Pa 
(np + a-2)(p'’P"“-'-a'’'’"''’-') 


( 2 . 19 ) 


The Bayes estimate of the variance of q can be obtained after 
obtaining the following expression. 

p 

E(e') = |e- 


N, 

QHp+a 


d0 


N, 




(np + a-3) 

a'P"'’"''! 


( np + a - 1) 


(np + a-3) 

Using these values in the following expression, variance can be 
calculate. 

v*(e) = E(e')-E-(e) 

. -E(0v)-(ef)' 


( np + a - 

'-P 

'a"’'’-"'] 

( np + a - 1 ) 2^-1^ “ 1’'-' ' 

,, .p- 

j 

(np + a- 

-3) p"'’"‘-'-a 

!!(> • ,1 - 1 

(np + a-2) P'"’"' ' ' 


.( 2 . 20 ) 


(-^ 0 ) 


The Bayes estimate of the survival function S(t) given by (2. 1 4) 


will be 


s-«) = J 


1-t 




N, 


0 np+a 


-de 




1 -TN 


} 


-p(n+r)-a+l __ ^~p(n+l)~a+l 


-p(n + l)-a + l 


Which has the following fomi: 
S*(t)=l- 


.(2.21) 


(,ip + a-l) [p'-N—i-a'"*'..-.'-' 


((n-l)p + a-l) 

The Bayes' estimate of the variance of S(t) can be expressed 
V*[S(t)] = E[S=(t)]-[S*(t)r 


Now, 






N, 

01. 


d0 


I’ 

= 1 + t-'’Nj0"'''’"^"-"d0-2l'’N| JO'""''‘^'’d0 




I - ps n »-2 i-a + l 


-p(n-f2)-a- 


- 2 l‘’N 




-p( n + ! ) - a -r 1 
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Which has the following fonn: 


E[S-(t)] = I + t 


,p (np + a - 1) [P 


( n+2)p+a- 1 _ n+2 ip+a-l 


((n-2)p + a - 1) a"'’P 
(np + a-1) [P*"' 


2p^(n+2)p+a~l _ i3r2)p+a-l 


!)p+a-l __ ^(n + Dp+a-i 


((n - l)p + a - 1) aPp' - P'W’"' 

From (2.21) and (2.22) putting the values in the expression of variance 
of S(t) we can obtain V*[S(t)] 

The conception rate function of the distribution given by (2. 1 5) is 


W(t) = 


pt 


p-i 


0P-t'’ 


The Bayes estimate of conception rate w(t) will be 


w (t 


pt 




N, 


fi 


0np. 


-d0=i6- 


np-a-p 


C) 


d0 






> 


p( n-f k -t-i I- 


dfi 


^,^p( U'-kt-! ipta-i 


k=() p( n -r k + 1) -r a - 1 


= pt'"'£t 


pk 


k=0 


(np + a - 1) 

’ j p( 11 k 1 1 T ; i - i _ ^ ^ p. I a r k • I p a - 1 

( n i“ k + 1 ) p + a — 1 < 

1 . k 1 ipiK n - k - 1 1 -a - i __ j:) t k H i ___ jy k • 1 s p^^ }i! 


Now, 


4 


E(w“(t)) = 


-np-a-2p 


^tV 




d0 


2t2,p-»y (k^ l)tPk 

^ (n + k+2)p + a-l 


k=0 


[p 


p(n + k + 2)+a~I __ ^p( n+k+2)p+a-l 


'] 




p(n + k + 2)+a-! n(k+2)p p{ n + k-f-2)p+a 




+ 11—1 


(2.24) 


• 


Using the values from (2.23) and (2.24) in the following expression 
the Bayes estimate of variance of conception rate can be obtained. 
V*[w(t) = E[w’(,t)]-[w’(t)]- 


The expected waiting time of first conception is obtained as 


E'(t) = j 


P .6-Ii-ci6 


Jp + ] 


Since 


E ( t ) = -!^0 

p + 1 

!' 

P + 1 

p + 1 (np + a-2) 


(2.25) 




And the variance for the same can be obtained by putting the 
value of 

E*(r)= 

Jp+2 

, a ^ 

Since 


E(r) = 


p + 2 


P + 2 


Nj0-np-n+2d0 


r.„., P (np.a-l)[a=P'’"--PV'"1 
(p + 2)(np + a-3) 

in variance format 


. V[E‘(t)} = E’(t-)-[E'(t)f 


2.3.2 A NEWLY DEFINED PRIOR DENSITY 

In this section a new continuous distribution is being used as 
prior density defined by probability density function of the underlying 
distribution is 


-g(e)^ 


he"'’ 

q” - IV’ ’ 


h<0<q; h,q>0 


(2.27) 



The Conditional likelihood function under the assumption of prior 
density (2. 1 6) provides the following posterior distribution 






H.p"0- 


q''-h” 


Which solves out to be in the followinn form: 




H 

j0-np.h-ld0 


where N: - — 


0-np+h-l 

h - np 

qh-np _ j^h-np 


( 2 . 28 ) 


The Bayes estimate of the parameter 0 can be obtained as 


0’ = 

J Qnp+h+1 


(h-np) 

_ Li ! p ->Q) 

(h-np+l) ’ 

The Bayes estimate of the variance of 0 can be obtained after 

obtaining the following expression. 

E(0C= fe-— f-d0 

' J /yip+li + l 

h ^ 

(h-np) 

“(h-np + 2)"q''-'’^'-!T-‘"’ 


.( 2 . 29 ) 



Using these values in the following expression, variance can be 

calculated. 

V*(0) = E(0')-E-(0) 


-E(0')-(0-)^ 

(h-np) q"-''P^ 2 _j^h-„p .2 
(h-np + 2) 


(h-np)q 


Ii-np+i __ I h-np+l 


(h-np + l)q"-"'’-h'’'"'’ 


.(2.31) 


The Bayes estimate of the survival function S(t) given by (2.14) 


will be 


S'(.)=j 


1- 


/tA'’ 


v0y 

q 


N. 


0np+h+l 


d0 


■tPNj0 


-(n+l)p-h-l 


d0 


h 




h-( n+I )p _ |.^h-(n+I)p 


h-(n+ l)p 
Which has the following form: 


S\t) = !-t 


„ (h-np) q'P'‘-"P-h 'p"-')P 


(h-(ii + l)p) q 


h~np |..h-np 


.(2.32) 


The Bayes estimate of the variance of S(t) can be expressed as 

v‘[S(tyi = E[s-(t)]-[s‘(t)]' 


Now, 
E[S^(t)] = | 






N, 


0np+h+ 


-d0 


= l-t"’Nj0‘'’"' 


2)p-h-l 


d0 


( 56 ) 



h-(n+2)p __ , h-(n+2)p M 

= i+t-PN,-2 9tpN^ r0-iN+i)p-i'-i(j0 

- h-(n + 2)p . -J 

Which has the following fonn: 

E[Sht)]= l + t-P- ,_ Jh-np) _ 

(h-(n + 2)p) q'’-'’P-h''-p (h_(n + i)p) 


■qh(ii + I)p __ |^h(ii + [)p 


(qi'-'>p_h''-''P) 

(2.33) 

" From (2.32) and (2.33) putting the values in the expression of 
variance of S(t) we can obtain V*[S(t)] 

The Bayes estimate of the conception rate given by (2.15) under 
the assumption of prior density (2.27) will be 


w7t)=f-^ 
’ J AP - 




0P _ p’ 0''p+i'+ 


pb’-'Nje 

h 

q 

pt'’''Nj0 


-np-Ii-p-l 


1- 


AP 


■p( n-1 )-h-I 


0 
{ t 


m 


SU) 


d0 


Solving the above expression we get; 


pt 


p-i 


1^ / I ! .»- t I \ 1 


h-{n+k --lip __ n-rk-^hp 


(h-(n + k + l)p) (q"-'-*'-h'-”b) 

Now, to obtain Bayes estimate of the variance of the same 


(2.34) 


E[w-(t)] = } 


li L 


pt 


P-1 




N, 


0iip+h+ 


-d0 


p-t-‘”^'’NjN,0' 


np-li-2p-l 


( t 




dO 
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f t A'”' 


k=0 




de 


Solving the above expression we get; 


= ph'<'’-'>£(k + l)tP 


(h-np) 


li“(n+k+2)p ___ n+k+2 )p 


k=0 


(h-(n + k + 2)p) 


h-np _ j^h-np 


.(2.35) 


Putting the values from (2.34) and (2.35) in following expression 
V * ( w(t))=E(w-(t)] - [ w* (t)] ^ 

The expected waiting time of first conception under the assumption 
of prior density (2.27) is obtained as : 


E-C) = l^e.^ 


d0 


vE(t) = -^e 

P + 1 


P+1 


N2j.e"’'’‘'’de 


p+i 

And 


h-np+l _ ] h-np+l 


h-np q -h 
(h-np + l) 


(2.36) 


ir(t") = | 


hM 


P_ 

:.p+2 


. 0 -.- 


N. 


0 '"’" 


11.-1 


d0 


p+1 


P + 2 


P + 2 


N J Q-np-h-ide 

h 

h-np 

h-np + 2“q"-'’'’-h”'""’ 


.(2.37) 
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Putting the values from (2.36) and (2.37) in following expression, 
we can get the expression for the Bayes estimate of variance of the 
expected waiting time of first conception. 

V*(E) = E’(t-)-[E‘(ttf 


2.3.3. INVERTED GAMMA PRIOR DENSITY 


The probability density function of inverted gamma distribution 


is; 


i(9) = 


/ 16 i 

7 V ) e 


v.i^.G > 0 


.(2.38) 


(0; rai+i) 

The applicability ofthis model in the Bayesian analysis is already 
being discussed by Chhikara and Folks (1971 ). The distribution reduces 
to exponential density for ).t = 0 

. The likelihood conditional under the assumption of prior density 
(2.3.8) provides the following posterior distribution. 


(y] 


l.p'’0 


HQ-np 






f V v“* 


0) 


r(ti + i) 


(e) 


ci-i e 


gJ 


r(d + l) 


d0 
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which solves out to be in the following form: 


G(0|tph hj = - 


Q-np-n-lg 1 9 


-(!] 

-np-p-lg Ui^0 




G(0|t„t, 


where N 3 - ' 


.(2.39) 


r(np + j.i) » 

The Bayes estimate of the parameter q can be obtained as; 


fi' 

0’ = j0.0-'''’'‘'"'e^”^N3d0 
0 

^Nj0-''P-.‘'-'.e^eJd0 

0 

which solves out to be as; 


0 *= (2.40) 

(np + j.i-1) 

The Bayes estimate of the variance of e is obtained after obtaining 
the following expression. 

E(O-) = j0h0'"''^“''e ‘‘’^N^dO 

0 

^ -1 i 

= N,j 0 ‘"'’-“''e"'’G 0 

0 

which solves out to be as; 

E(e=) = (2.41) 


(np4-p-2)(np + p-l) 
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Putting these values in the expression variance can be calculated 


as 


V'(0) = E(e-)-(9*)' 


(np + |.i-2)(np + ia-l) 
V = 


(np + p-1) 


.(2.42) 


(np + p-2)(np + p-l)- 
The Bayes estimate of survival function under the assumption of 
prior density (2.25) is obtained as 


s'(t) = l 


t Y 

V0V 


j0-np-MMg 






‘ '■ ‘3 „(n+l)p-.i. 


Which has the following form: 
, ,, ,, rf(n + ! )p + p] 


.(2.43) 


r{np + LL) 

The Bayes estimate of variance of survival function can be 

expressed as 

V*[S(t)]=E[S(t)]-[S‘{t)!' 


where 

■r> 

E[S-(t)] = | 




e; 






0 

- ':^^p ^,p r[(n + l)p + d] 

I (np + j.i) r(np + d) 

The Bayes estimate of conception rate given by (2. 14) under the 

assumption of the prior density (2.27) is obtained as 


.(2.44) 


w’(t) = J 


„fP-' 


Pt’ 


e^-tf 


e-'p-p-'e ^«^N,de 


. pk 


= pt^-'N, f - y - do 

'o (0y^(0) 


k-n 

Solving the above expression we get; 


\v 


*(t) = Dt‘'^'N y npCn + k + h + p] 

■ r[(np + 5) 


and also, 



k vO 

*! 


_ fl 

r pt-' 

-J 

0 1 

0' -f' 

■‘J 

'N. 1 



0 -..p-M-ie i<’y,d0 


p:^2,pM,^,_Jj) P(n. 2 ,-d-le" + 


( t 


d0 


k^O 


J 


.(2.45) 


Solving the above expression we get; 

,(k+ 2 )|. r[p(n + k+2) + p| 


E*[w-(t)]-ph'"’-''’X(k + lKt'’y 


k=0 


r[(np + 6) 


..(2.46) 
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Now, the Bayes estimate of the variance of conception rate can be 
obtained by putting values from (2.45) and (2.46) in the following 
expression. 

V*[w(t)] = E(w-\t)]-[w*(t)]^ 

The expected waiting time off first conception under the 
assumption of prior density (2.25) is obtained as; 

=0 f-V f n ^ 

E*(t)= f— vE(t) = -L .0 

ip + 1 ■ I P + 1 )' 

f-1 

p+ii 

_ P V 

p + 1 np + p - 1 


(2.47) 


And 

y.' ( '■' ) 

f_P_.0he-"P-p-'e2«lN,d0 

f.P + - 


r; 



“^N,d9 


V l-(t') = 



A 

,/ 


p + 2 ( np + p - 1 )( np + p - 2 ) 


( 2 . 48 ) 
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Now putting the values of (2.47) and (2.48) in the following 

expression, we can obtain the variance for the same as; 

. V*(E) = E*(t^)-[E‘(t)]^ 

or 1“ 

p P ^ 

p = 2 (np + (.i-l)(np + g-2) [p + l (iip + E^l) 


pv- 


(np^p-1) 


1 


(p + 2)(np + p-l) 


P 

[(p + l)(np + p-l)]- 


(2.49) 


CHAPTER - 3 

STUDY OF BURR TYPE III 
DISTRIBUTION USING TYPE I 
CENSORED DATA 

3.1 Introduction 

Generally there are two type of censoring Type I and Type II 
censoring. In Type I censoring, the number of failure is random whereas 
the length of the experiment time T is fixed. In Type II censoring, the 
number of failures r, is fixed while the length of experiment time T is 
random variable. 

This study considers Burr Type III distribution, introduced by 
Burr ( 1 942) as a versatile unimodel distribution. The Burr Type III family 
is much richer than Burr Type XII as is shown by the moment ratio 
diagrams constructed by Rodriguez (1977, 1982) Burr Type 111 
distribution includes exponential, Weibull and Log-logistic distribution 
as special cases. The curve shape characteristics for Pearson Type I 
(Beta), II, V, VII, IX and XII families (see Rodriguez (1977, 1982) is 
also present in Burr Type III distribution. As the model based on this 
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distribution is very rich and simplest functionally therefore it is very 
useful and most attractive for statistical modelling. Wingo (1993) has 
described method for fitting the Burr type XII distribution to complete 
life test data by maximum likelihood method. Recently, Wingo (1993) 
consider Burr Type XII distribution and discussed Maximum Likelihood 
method for fitting the Burr Type XII distribution to progressively censored 
life test data. 

The present investigation extends the Wingo (1993) work by 
considering Burr Type III distribution, which is much richer than Burr 
XII. It also extends the work of Wingo ( 1 983) by developing mathematical 
and computational methodology for fitting the Burr Type II distribution 
to type I progressively censored life test data. Mathematical expressions 
are obtained for finding the asymptotic variances and covarience of the 
MLEs of the parameter of the distribution. The existence and uniqueness 
of the MLEs for arbitrary sample data are also discussed. 

3.2 Models and Assumption 

Let {Xi>0J = U2 ,fi}be the independent lifetimes random 

variable having a two-parameter Burr Type III distribution with the 
cumulative distribution function F(x) = [l + x'n'‘'; x^c,k>0 (3.1) 
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Where c and k are the shape parameters. Thus the probability 

d 

density function is — F(z) = + 

ax 

^ (1 + x'^)'^"' (^-2) 

Putting k=l , Burr III reduced to log-logistic with probability density 

function is given by f(x) = — — — yf (3.3) 

The Hazard rate function of Burr III distribution becomes 

: (3-4) 

JcoT^ 
l-(l + x'^)'' ■ " 

, , , kcx'''"'"'’ 

(3.5) 

and the quartile function of the Burr type III distribution is 

Q(P) 

p= jf(x)dx 

0 

F = (Q(P)) = P 

[i+(Q(prr'=p 
i+(Q(pr = p"‘'^ 

(Q(pr=p'''^-i 

(e(P))-'=y«-' 
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1 


(Q(P)“^ 



J/k 


Q(P) = 


p./kc 




Suppose that the total number of items (N) are placed life testing 
experiment and let n be the number of known exact failure. Further suppose 

that the censoring occurs j = 1,2,3 ,m. If, at j”'*’ stage of censored the 

randomly selected items r^ are removed or censored from the test. Thus 
we can write. 

m 

N = n+Zrj (3.6) 

j=l 

and the likelihood function for m-stage Type I progressively 
censoring is given by 

n m 

[J=xnf{x,)n[i-F(T,)]'' (3.7) 

n=l j=l 

for some constant x ■ 


3.3. Maximum Likelihood Estimation 

The log - likelihood function (LLE) based on the I-state Type-1 
progressively censored sample, can be written as 
L - In ZJ = n\ogk + n\ogc-{c+ I).Llog.v,- - {k + 1 ) Zlog(l + A' ‘ ) + S log[l -(1 + Tj^)-k} 

/=! /=1 7 =! ■ 

(3.8) 
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Differentiating partially the above equation with respect to c and k 

respectively, we get the following likelihood equations 

5L_n " , " xr'InXi (l + T:D'^'"'’VlogT, 

dc c i=i i=i(l + xr‘') j,/ + (j-y) 


5L n 


" kr.('l + T''r' ' 

--Zlog(l + x-)-Z -i" ^ ^ =0 

5k k i=i H [1-(1 + T:-^)'''] 


(3.10) 


The MLE's gand f of c and k respectively are the solutions of 
the above equations (3.9) and (3.10). These equation can not be solved 
explicitly. But they may be solved by using numerical method in order to 
determine the values of c and k. 


3.4 Asymptotic Variances and Covariance's 
ofMLEs 

The observed information matrix F, is obtained by taking negative 
of second partial derivative of log likelihood function L with restart to c 
and k. Thus 

-S-InL. T 
6c§k J 

-a-InLj 1 

ak’ 


n 

Z 

i=l 

-a- logLi 
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"-a^ logM" 
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1=1 
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F = 11 


f f 

^00’ •‘-0: 
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Where fgh = 


5-L 


50g 50n ’ 


g,h = 0,l 


(3.11) 


3.5 Existence and Uniqueness of MLE^S 

Case I : Suppose k is known and fixed, we have to show that the 




estimator c exists and is unique. For uniqueness, it s to show that the 


function defined in (3.8) is strictly concave in c. Since 


is negative, it 


implies that the function is strictly concave in. Hence c is. For existence 


of c, we have to show that the equation 


zero, finite zero on the positive real line, provided that Xi * 1 and / or 

Xj?il for some i=l, 2 ...., n and j = 1, 2 , m. 

Let us define eq. (3.9) as a function of c only and index sets as 


‘ r.THnT; 



As the observed asymptotic variance-covariance matrix of 
estimators £ and f is the inverse of the observed information matrix F, 


therefore 


var(c) cov(c,^) 
co\{k,c) var(^) 


After some algebraic simplification, we obtain 




{k + 1)(1 + ^ + 1)(-1 + j-cY2(k^\) _ + j’-c) 


■ 2 (^ + 1 ) 


{k + 1)(1 + + (1 + jjcyKM) 


After some algebraic simplification and direct calculation under 
the above non-empty index sets, as c -> 0 and as c -> oo , equation (3.15) 
yields respectively {j)(0) = oo 

(j)(oo) = k E log Xj - E log Xj - k E rJnTi H 1 V"! 

i = M i = P ^ J w 


The right hand side of (3. 1 7) will remain negative, as long as either 
of the index sets M, P, Q or R. Therefore under these conditions on the 
sample data we may say that (j) (0)<oo and (j) (oo)<0, there exist at least 
one positive real line. 

Differentiation and straight forward algebraic simplification show 

d(j)(c) 

that (i.e. ^ (c) is monotonic decreasing function in c), implying 

this root c is unique. Therefore we conclude that (i)(c)=0 has a single 
finite zero on the positive real line. Hence the MLE's is always exits and 
is unique for known and fixed k. 

Case II : Suppose c is known and fixed, the function in (3.8) is 
strictly concave in k, because eq. (3.14) is negative. It follows from the 
strict concavity of (3.8) that k is unique now we have to show that k 
always exit possesses a single, finite zero on the real line provided that 

Xj?i 1 and / or T. 1 for some i = 1, 2, , n and j = 1, 2, ...., m. 

We define 

n " " . ' . ' ri[(l+r)Mog(l+r)-cTffnT,] 

4^(k) = -+cEiogXj - E log(l + xf ) - Er|(l +1)' ) + E^ n+Tf 
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After some algebraic manipulation and direct calculation of ecj. 
(3.18) over non empty index set defined in eq. (3 . 16) as k-^ 0 and k oo 
we get (°°) = -c i log(Xi) - 1 log(l + xi-^ ) - I (1 + t; ) we see that right- 
hand side of the eq.(3.l0) is negative if and only if index set Q or R is 


non-empty (i.e. (!)(oo) < 0) . thus <t)(0) < co and (t)(oo) < 0 , implies that there 


exist at least one real root of ^(k) = ^ ^ • one can easily show that of 
< 0 . i.e. is monotonic decreasing in k. 

It follows that this root 's unique Therefore we conclude that 
the eq. T(k) = 0 has a single, finite zero on the positive real line Hence 
when c is known and fixed, the MLE k always exists and is unique. 


3.6 Conclusion 

This chapter develops mathematical methodology for estimating 
the parameters involved in Burr Type HI distribution from Type 1 
progressively censored life test data by maximum likelihhod method 
Asymptotic variances and covanence's of MLEs are obtain. We obsci ve 
that the MLEs of the parameter exist, unique and finite for di bin ary 
sample data. We also observed that the result of this chapter be applied 
to type I progressively censored sample data arising in a life test 
in clinical trials and fields where Burr Type III distribution is 


experiment in ( 
used. 
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CHAPTER - 4 


A STUDY OF GUMBEL'S BIVARIATE 
EXPONENTIAL DISTRIBUTION FOR 
LIFE TESTING ANALYSIS 

41. Introduction 

The most widely used and accepted model for life testing is the 
exponential model. The density function is given by: 

/(6,t) = 0exp(-0.O; 0,t>O 

The two important properties of this model viz. Forgetfulness 
and constant hazard rate make this model even more popular in life 
testing analysis. 

The constant hazard rate means that the exponential distribution is 
an appropriate model for the life time of an item when there is no ageing 
or wear out. It is mostly applicable to electronic components since they 
fail due to chance and not due to wear out or ageing. This model is 
equally popular in the similar branches of studies like fertility and survival 
analysis. 
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The utility of it, leads a number of workers to develop the bivariate 
(multivariate) extension of this distribution. There are a number of bivariate 
exponential distributions. 

In the present chapter an attempt has been made to characterize 
the concomitants of the order statistics of Gumbel's exponential bivariate 
distribution. 

Let (x^y,), i = 1, 2, ...., n be n pairs of independent random 
variables from Gumbel's exponential bivariate distribution (which may 
be consider as a particular case of Gumbel's Weibull bivariate distribution) 
with the distribution function. 

F(x,y) = l-exp[-ax]-exp[by] + exp[-ax- by -cxy]; 

(x,y>0),(a,b>0),(0<c<ab) (4.1) 

and the corresponding probability density function can be obtained 


as; 


r(x,y) = ^{^[F(x,y)l 


dy 


{a.exp(-ax)-(a4-cy).exp[-ax- by -cx>']} 


= {(a + cy )( b + cx) exp[-ax - by - cxy 1 - c exp[-ax - by - c.xy ] } 
which finally yields; 

f(x,y) = [b + cx)(a + cy) - c]exp(-ax - by - cxy); 

(x,y>0),(a,b>0),(0<c<ab) ...... (4.2) 
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4.2 Marginal Probability Density Functions 

In the present section the marginal probability density functions 
and distribution functions of X and Y are obtained. 

The marginal probability density function of X can be obtained 


as; 

00 CO 

g(x)= 

-oO 0 


= I [(b + cx)(a + cy) - c] exp(-ax - by - cyx)dy 

0 


oc 

= j [ab + bey + acx + c^xy - c] Qxp[-yib + cx )]dy 

0 


= e 


-a.\ 


I cy=-'ib + + {ab 

,0 


+ acx-c)\e-y^^^^-^dy 
0 


Since j exp(-p0t''dt = p-‘'’"'T(v + 1) 

0 

-as[ c(b + cx) (ab + acx-c) 
|(b + cx)-^ (b + cx) 

c + ab-c + acxl 
(b + cx) J 

I (b + cx) 



: ae 


which finally yields 
g(x) = a exp(-ax) ; x>0 
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( 4 . 3 ) 


And the marginal distribution function of X is; 


G(x)= jg(x)dx 

0 


= |a exp(-ax)dx 
0 

G(x) = l-exp(-ax) ; x> 0 (4-4) 

Similarly, the marginal probability density function of Y will be; 

CO CO 

h(y) = j f(x,,y)dx = j f(x,y)dx 

-00 0 


CO 

= j[(b + cx)(a + cy) - c]exp(-ax - by - cxy)dx 
0 . ^ 

00 

= j [ab + bey + acx + c'xy - c] exp[-x(b + cx)]dx 

0 


CO 

= e’*’’ j {x{ac + c V) + {ab + bey - c)] exp[-A-(fl + cy)]dx 

0 

_^-bv[ b(a + cy) | 

1 (a + cy) J 

which finnaly yields 

h(y) = b exp(-by) ; y > 0 (4.5) 

The mean of the random variable Y having the probability density 

function as defined in (4.5) will be; 

00 

0 
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00 

= Jybexp(-by)dy 

0 

G(Y)=i 

Also; 

00 

E(Y2) = \y^Ky)dy 

0 

00 i 

= ^y^btx^{-by)dy 
0 

E(Y^) = 

And the marginal distribution function of Y is; 

V’ 

Y[{y)^]Ky)dy 

0 

y 

= Jbexp(-by)dy 

0 

H(y) = l-exp(-by) .; y > 0 


(4.6) 


(4.7) 


(4.8) 


4.3 Conditional Probability Density 
Functions 

Now, in the present section the conditional probability density 
functions of Y for given X and of X for given Y are obtained. 

The conditional probability density function of Y for given X can 
be obtained as follow; 
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(4.9) 




f{x,y) 

gix) 


/( 3 ;|x) = [(^ + + cy) - g] exp[->’(6 + cx)] 

a 


Similarly, the conditional probability density function of X for 


given Y will as follows; 


f{x\y) 


f{x,y) 

h{x) 


[(b + cx)(a + cy) - c] exp(-ax - by - cxy ) 
bexp(-by) 


f(Ay) 


{{b + cx){a + cy) - c] exp(-.x(a + cy)] 
- 


(4.10) 


4.4 Probability Density Function of Order 
Statistics From Gurabel's Exponential 
Bivariate Distribution 

In the present section, the probability density function of the r'*' 
ordered statistics and the joint probability density function of r and s 
ordered statistics from Gumbel's exponential bivariate distribution are 

obtained. 

For Gumbel's exponential bivariate distribution the probability 

density function of the f'' order statistics X_..^ is, 

4(x) = C„[F(x)r'[l-F(x)rf(x); 0<x<oo 
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Where = 


n! 


(r-l)!(n-r)! 

= frn = - exp(-ax)h“’ exp[-a(n - r + l)x] ....(4. 1 1 ) 

In particular for r = 1, i.e. the probability density function of the 
first order statistics is; 

frjt (^) ~ exp(-«(2x) (4.12) 

For the distribution with probability density function (4, 3), the 

joint distribution of two order statistics and s'*’ is as follows; 

= Q,JF(x,)r'[F(x2) - F(xOr'--ai - F(x2)rY 

where = • 


n! 


(r-l)(s-r-l)!(n-s)! 

= Q,.:«[l-exp(-ax,)r' 

{[1 - exp(-ax 2 )] - [1 - exp(-ax, )] 

[expC-ax, )l''‘'[a exp(-aX| )][a exp(-ax, )] 

= C, , „[1 -exp(-ax,)]^'' {exp(-aX|) -expt-ax,)]}^'' 
.[exp(-ax 2 )]'’‘'a' exp[-a(Xi +X 2 )] 


.s- r - 1 


1=0 i-0 


'r-f 

's-r-l 

i 

j 


exp(-ia + sa - ra - a - ja)x,)] 


[exp(-jaxj )]exp[-(na - sa)x Jexp[-a(X| + x, )] 

which finally yields 


r-1 s-r-1 


i=0 j=0 


'r-f 

's-r-f 

i 

j 


exp(-(ia + sa- ra- ja)x,)lexp(-na-sa + ja + a)x,] (4. 13) 

exp(-(i + s- r - j)ax, ) exp- (n - s+ j + fiax, 
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. Yr-lYs-r-H 

C„.a=SE(-')'*' 

i-0 j=0 


i A J 


4.5 Probability Density Function of 
Concomitants 

Now, the probability density fiinction of the first order concomitant 
(i.e. r=l) of the order statistics is 


0 


I [(b + cx)(a + cy) - clexpl-y(bH^ na exp(-nax)dx 

i a 

A 


= n j [(b + cx)(a + cy) - c]exp[-y(b + cx)] exp(-nax)dx 

0 


= nexp(-by) j [cx(a + cy) + ab - c + bey] exp[-x(na + cy)]dx 

0 

Now using the following relation, 

00 

jexp(-pt)p''dt = p'‘'’"‘'r(v + l) Re(v)>-l,Re(p)>0 
0 


= nexp(-by) 


c(a + cy) ^ (ab-c + bey 
(cy + na)‘ (cy + na) 


After simplification. 


gti:,.)(y) = exp(-by) 


nb- 


abn{n-l) an(n- 1) 


^ na ^ z' 


Y + - 

V' c ; 


na 

y + - , 

V C J 


V 


;y>0 


(4.14) 


Similarly, can be obtained as; 
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g(i:n)(>^) = exp(-by) 


nb- 


y + 


nbY 

c J 




■> 

bn(n-l) 



( 

T J 

c 

x + — 



( c ) 



;y>0 


.( 4 . 15 ) 


Now, to prove that g(|.„)(y) is a probability density function, we 

have to prove the following identity; 

00 

jg(l:n)(y)dX^l 

0 

Taking Left Hand Side (LHS) 


CO 00 

jg(i:„)(y)dx-|exp(-by)nb- 


abn(n - 1) 


na^ 

c 

y + — 


( c J 


c y + - 


na 


dx 


00 

= nj bexp(-by)dx-| 


exp(-by)dy 


( aa) 

1 ( na^ 

c y + — 

c y + — 

V c j 

V c J 


Let 


na ) 
V + — 

V c J 


= t 


then above expression becomes 


: ' 


Using the relation 

00 

je"“u'‘du = E|(ax) 
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We get 

OD 

Jg(i:n)(y)dy = ri-(n-l) = l 

X 

Now, the probability density function of the concomitant can 
be obtained by using the following formula as; 




■n+r~\ 


f i-lYn\ 


/=/i=r+l 




gihn(y) 


where 


g(i:i)(y) = exp(-by> 


Thus, 


abi(i-l) ai(i-I) ^ 


la 

y + — I c| 
c 


f ia 

y + - 

V c 


.( 4 . 16 ) 


«(r,,|(.V)= 

/=//~r+l 


/-/?+/*-! 


f 




exp(-6y) 


ib- 


abi{i-\) ai{i-\) 


Y + 

\ c 


ia ^ f 


la 


V c 


Similarly we can obtain, 


g(r:n)= 


i-n+r-l 


i = II - r + 1 


vn-i'AV 


exp{-by) 


ib- 


abi(i - 1) 


ai(i-l) 

( ia^ 


( ia Y 

c x + 

1 cj 

x + 

V c J 


l c ) 
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Now, the moment about origin of will be, 






i=n-r+l 


Vn-ryw 




Where 


, ^.r)k + l) fia 


Hence, 


vcy 


(i-i)r(k+i)expfi5^Yr(i-k,i^ 
V c ;v c 


St-'r 


n - 1 Ynt 


vn-ryv; 


T(k + 1 ) fiaW; 




— J (i - l)r(k + l)exp|^— r|^l-k ^ 


iab' 


(4.20) 

Similarly, the k* moment about origin of X . can be obtained as 


S(-') 


i-n+r-1 




1 = n -- r + 1 


lYn) 


V^-^AV 


■ r(k + i) (\h)\. 


c j 


(i-l)r(k + l)exp 




r 


C / 


iab 


r(i-k, 

V c 




(4.21) 


Now, in particular for k=l the mean of will be; 


E(Y,„|) = |X„„,= Z(-l) 


IL i-n+r-l r i -1 YnY 


i=n-r+l 


In-rAi) 


1 la,. 

(i-l}exp 

I b c 


' iab ^ 


fn iab)l 


r 

o.~ 

c j 


c jj 


.(4.22) 


The variance of Y(r:n) can be obtained as 

= th'(r;n) ” Xl-htrm) 

Where 




C(- 

i=n~r+l 


^ i_l\YnY 


(n-rAi; 


2 J ia 




KcJ 


(i- l)exp 


V c ; C 


.(4.23) 
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Similarly, the mean of X(r:n) will be; 




/■-/i+r-l 


i=n-r+\ 


fi-\Yn\ 


yn-rj^ij 


i ib 


a c 


{i - 1) exp 


V c j 


iab 


0 , 


.(4.24) 


Teh variance of X, , can be obtained as 

(r:n) 

u 

Where 




i-n-r+1 


f i-i A 


^n-r^ 




VV 


-21 

a- V c ; 




^ iab^ 




iab 


.(4.25) 


It may be noted that; 


I^U)=n^^| = nE(Y) 


.(4.26) 


and 


f 2 ] 




r = l 


.(4.27) 


VD J 


4.7 Joint Distribution of Two Concomitants and 

Now, for the Gumbel's bivariate exponential distribution with 
probability density function (4.2) the joint distribution of two concomitants 


of order statistics Y . and Y , can be obtained by using the following 

(r.n) (s.n) 


relation; 


fy{nnXy{s\n )(y\y2) = g(r,s-.n)iy\^y2) • 

C0X2 

Sir,s:n)iyiy2) = j J/(3^ll^l )/(>'2l^2)/r,^;„ (-^1 Al 


0 0 


Where f ^._^(x,X 2 ) is the joint probability density function ofr^*" and 
s'’’ order statistics of X. 

00 X-. 

=jli(b + cx, )(a + cy I ) - c]exp[-y , (b + ex, )] [(b + cx, )(a + cyj ) - cjexpl-y^ (b + cx, )] 


0 0 


r-l s-r-l 




i=0 j=0 


f r-lYs-r- 1 ^ 


exp[-{ia + sa - ra - ja)x, ] exp[-{na - sa + ja +a)x2 JdXjdXj 


(b + cXjjCa + cy,) - c] exp[-(cy2 + na - sa + ja + a)x2 ){I(x, ,y, }dx, | 

(4.28j 


Where, 


■'2 

I(x,y,) = |[(b + CX, )(a + cy , ) - c]exp[-cy , + ia + sa - ra - ra - ja)x, ]dx, 


Integrating by parts and simplifying the above expression, we get; 
- J [b(a + cy , ) - c]--^- (exp[-k(y , )x, ]-\\ 

k(y,) J 


with k(y,) = cy, + ia + sa-ra-Ja 

Putting the vlaue of ICx^, y,) in equation (4.28), we get, 

r-i s^i /r-lYs-r-l) 

=:C,,„exp[-b(y,+y2)]i; 

i=0 j=0 
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I i A J ) 


r-1 s-r-1 


[(b + cxj )(a + cy j )]X E (“ l)(i + i)| 

i=0 j=0 


rr-lYs-r-n 


A j J 


c(a + cy|) 


. I X2exp[-k(y,)x,]--— exp[-k(y,)x,]-— 
k(yi) i k(y,) k(y,)^ 

|[b(a + cy , ) - {exp[-k(y , )X2 ] - 1}| 


Hence we get 

. / r-lY s-r-H 

g(r,s:n)(yiy2) = C,,„exp(-b(y,y2))E 

i=0 J=0 


J A J ; 


(a + cy,)(a + cy2)j 2 


1 


1 


+b'c' 


' k(y,) [' [k(y|)kYy2) k(y,)k(y2)^ k(y,)(k,(y|) + ktys) 

k(y,) 

k(y2)(k(y,) + k(y,) 


+ - 


cYa + (y,) 

f 1 , ' 

_ ' 1 

k(y,) 

Yk(y,) + k(y 2 ))- k(y,)(k(y,) + k(y,) kCy^kCy,); 

bc((a + c’^,) 

k(y,) 1 c-{a + (y,) 

1 ' - ' 1 

k(Y,) 

k(y:)(k(y,) + k(y,) k(y,) 

Uk{yi) + k(y:)' k-(y,)J 

bc(a + (y2) 

k(y,) 

1 

C' 

k(y,) 

’k(y,)(k(y,) + k(y,) k(y 2 )(k(y,) + k(y,)) 


with k(y j) = CYi + ia + sa - ra - ja and 
kCy^) = cy, + na - sa + ja + a 
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CHAPTER - 5 


BAYESIAN ESTIMATION OF 
SYSTEM RELIABILITY 

5.0 INTRODUCTION 

In the context of modern technology and its future developments 
system reliability assessment is of paramount importance. The objectives 
of a system engineer include the study of various reliability characteristics 
of the system for improving upon its performance. These characteristics 
can be better analyses if apart from the observed life data the 
environmental experience of the variations in the reliability parameters of 
the system is also taken into consideration. Bayesian analysis combines 
the prior belief about the life time parameters with experimental data. 
The studies by Balagurusamy ( 1 984), Sinha ( 1 986), Mann et a! ( 1 974) & 
Kapur & Lambason (1977) used classical and Bayesian inference 
technique in case of some important life time. In many practical situations 
the operational data on the complete system are either limited or non- 
existent. In such cases, one makes use of the life time data recorded on 
its components. The studies by Kaplan et a! (1989). A apostolicism 


(1990) Martz and Waller (1982) reviewed the relevant literature on 
probability safely assessment (PSA) of a system and laid out a language 

conceptual framework andmethodologyfordealingwithsuchsituations. 

Kaplan et al (1989) analysed various probability curves useful for 
expressing the degree of confidence about the complete system reliability 
and the way in which Bayes theorem updates prior probability cuives. 
Highlighting the importance of posterior estimates foi leliability 
characteristics of various system networks. Shanna and Bhutan! (1992a) 
presented Bayesian reliability analysis using the truncated failure 

infoiination. 

For improving the performance of a system maintenance is used 
by reliability practitioners maintainability of a maintained system is as 

importantasreliabilityitselfforanengineer.Theconceptofavai, ability 

combinesbothreliabiUtyand,nainability,lnferencesonsystemava..ab,lity 

in detail in the literature both from classical and Bayes 


have been studied in ' 


view 


point. Most of these 


studies used simple probabilisitcs reasoning 


for analysing availability of various systems. Further, a system may vai y 

ry complex one. One approach for analysing such a 

of convenient sizes, each 

reliability and availability of susbsystems 


from simple to a vaiy 

system is to decompose it into subsystems 


repr 


resenting a specific function 
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are then estimated and combined to determine the reliability of the entire 
system using certain probability laws.- In the light of above discussion 
the present Chapter develop Bayes estimators of the system reliability 
and available in case of some important and basic system configuration. 

5.1. STATISTICAL ASSUMPTIONS 

For performing the statistical analysis on the system reliability 
characteristics 

(1) It is assumed that the failure time distribution for each 
component is negative exponential with p.d.f. 

f(x;0) = ee-'’^x,0>O (5-1) 

Therefore 0 and ^ are the failure rate and mean time to system 
failure (MTSF) respectively. 

(2) It is assumed that the repair time distribution for each component 
is also negative exponential with p.d.f 

. f{y;0') = 0'e^^’fy,0>o (5-2) 

Here and ^ stand for repair rate and mean time to system 
repair (MTTR) respectively. 
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(3) The reliability function for the system or compone 
may be defined as 

(5.3) 

R(t) = P[x>t] = e-^t >0 

Where t is called the mission time 


•1 u'litv is the probability that 
(4) The system or component availability ii> t' 

the system is operative in the long run and may be e>^P 

MTSF 


Availability = 


MTSF + MTTR 


0 


1 

e ef 


e+e' 


.(5.4) 


5.2. RELIABILITY AND AVAILABILI 
FOR SOME BASIC SYSTEM 
CONFIGURATION S 

• lowMoeen discussed earlier. 

Various basic system configurations ha 

ran expi'oss reliability 

It is notable that using basic probability laws one 

, u-iitates of its components, 
of the whole system in terms of the rehabi 

, • expi'ossions tor system 

Almost on the same lines, one can obta 

it has been assumed 

availability in terms of component availability. ’ ^ ^ 

■ f -ndeoendent and identical 
that each system configuration consists o i 
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components. Denoting the reliability and availability functions for series- 
parallel, parallel-series and a k-out-of-m system by 
l^sii,(^h'^sp,^ps( 0 .Ap 5 _,R^^(t) and respectively, some well known 
results are stated below. 


5.2.1 SERIES-PARALLEL SYSTEM 

This system consists of n-subsystems arranged in series and each 
subsystem has m-parallel components. 

For this system 

and 

Considering 


.(5.5) 


(5.6) 


a=0 


yaj 


Where 




yaj 


0, a > n 


the equations (5.5) and (5.6) may also be rewritten 


as 


a:=0 \U() 

and 

n im 

a„=ES<-') 

(i =--0 


a+^ 



^ ma'’ 

U/ 

1 P J 


|[R(t)f,t>0 


.(5.7) 


a+li 


fn] 

'"ma^ 




A' 


.(5.8) 
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5.2.2 PARALLEL SERIES SYSTEM 

In this configuration there are n-parallel subsystems each having 
m-component amanged in series. Here, 

R (t) = l-[l-{R(t)}'"]'’ .......(5.9) 


or 


-1 


a = l 


and 








[R(t)]'"“,t>0 


.(5.10) 


.(5.11) 


or 


A,.=E(-ir' 






yaj 


A" 


.(5.12) 



5.2.3 K-OUT-OF-M SYSTEM 


With simple probabilij'tic reasoning, the reliabililt for this system 


is given by 


[R(t)f[l-R(t)] 


5.3 DEVELOPMENT OF POSTERIOR 


DISTRIBUTIONS 


Now, since, the failure time distribution for the components is 
negative exponential with failure rate p. therefore, fora total testing time 
T, the probability of observing r failures will be given by the poison 
distribution with a p.m.f. 

(5.16) 


= 0 , 1 , 2 . 


r! 


Similarly, the number of component repair(s) perfonned in the me 
(0, 1 ') will also follow the poison distribution with p.m.f 

(5.17) 


. P(si;V) = .^-^-^^;s = 0,1,2... 


Now, the posterior distribution of p in respect of the prior in 
(5.5) given that r failures have been observed in the time interval (0, T) 
can be obtained as 

P(i1l-0g(p) 


n(|.i|r): 


j P(rip)g(.u)dp 




r{r + a) 


;|.L,(T + l),(r + a)>0 


.{5.18) 


Which is Gamma density with parameters (T+ 1 ) and (r+a). 
Similarly, the posterior distribution of /, in respect of the prior 


(5.6) can be put as 


n(>tis) 


P(s|^)g(^) 

GO 

jp(slX)g(A,)dX 

0 


r(s + b) 


;L(T + l),(s + b)>0 


(5.19) 


The distribution in (5. 1 9) is also a gamma density with parameters 
(T-Pl) and (s+b). 

Now, on using the respective posterior distributions of \x and x 

X 

in (5.18) and (5.19) and using the transformation ^ = ^7^ in (5.8) the 
posterior distribution of A given r and s comes out to be a beta 
distribution with p.d.f 


r(Air,s) = 


1 


B(s + b,r + a) 


(s+b).(i- + a)>0 (5.20) 


5.4. BAYESIAN ANALYSIS 

5.4. 1 . BAYES ESTIMATES OF RELIABILITY FUNCTIONS 
for SERIES-PARALLEL, PARALLEL-SERIES 
AND K-OUT-OF-M SYSTEMS 

For series-parallel, parallel-series using (5.5), (5.7) and (5. 1 9), the 

Bsyes point estimate of the reliability function for a series-parallel system 

can be defined as; 

R;,(l) = E[R,.(t)lr] 

= R,„(t)7t(t)Tx(nlr)di.t 
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a=0 p=0 

ni ma 

a=:Q fi=0 




yCX; 


^ ma^ 

vaA P ; 


(T+lp” 
r(r + a) g 




dj-i 



1 

f 

TR 


1 + - 


(5.21) 

A T + U 

and on using (5 . 5 ), ( 1 5 . 1 0) and (5.19), the Bayes point estimate of 
the reliability functionAora parallel-series system becomes: 


Rps (0 = E[Rp,(t)|r 

^ "At+1)'"" 


a=0 \^J 


£(_i)- 

a -■ 1 


r(r + a) 




~1 -(T+imt+l) 

^ FAlu 




fn) 

i 

v«J 

, mat 

1 + — 

T+i; 


Similarly, on using (5.5), (5.14) and (5. 1 9), the Bayes point estimate 
ofthe reliability function for a k-out-of-m system can be put as: 


III--' 


ll=t. v = 


m V I'll - u 


u j\ V 


V J 


(T+l)'- 


r(r + a) 




:i-l ut'VHI 

e 


dp 




u=r v=U 


VI 1 


{ V 


1 + 


(u + v)t 


(5.23) 
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5.4.2. BAYES ESTIMATES OF AVAILABILITY OF A 
SERIES-PARALLEL, PARALLEL-SERIES AND 
K-OUT-OF-M SYSTEMS 


of a series-parallel system can be obtained as; 
< = E[Ajr,s]; 

= |A3pf(alr,s)dA 


hhr ' UJU jB(s+b,r + a)J 

B(s+b,r + a) 


P J B(s+b,r + a) 


■ Similarly, on using (5.6), (5.12), (5.20), the Bayes estimate of 
availability of a pai'allel-series system can be found as: 


A:. = E[AJr,s] 




J A dA) 


B(s + b.r + a) 


B(s+b+ma,r + a) 


(5.25) 


[aj B^+b-r+a) 

Also, on using ( 1 5.6), (5. 1 5) and (5.20), the Bayes point estimate 
of availability of a k-out-of-m system becomes 


A';., = E[A,Jr.s] 


u=kV>-‘y 


u jB(s+b,r + a) 


j /^^s+b+u-l^j _ A)''^“*'’’''‘'^ 'dA 


"■y m ) B( s + b + u, r + a rm u) 
Z u B(sTZmt-a) 


.(5.26) 
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5.4.3 An Example 

For the purpose of studying the behaviour ot posterioi oi Bayes 
estimates of reliability and availability functions of series-parallel, paiallel- 

series and a k-out-of-m systems, we consider an example with m n 

k=3 and mission time t=30 units. Further also let i a and s 
recorded failure and repair information in the time interval (0, 100). Using 
(5.21), (5.22). (5.23). 5.24), (5.25), (5.25) and (5.26). the Bayes estimates 
of reliability and availability for varying mean failure rate 'a' and fixed 
mean repair rate 'b' have been tabulated in table-I and II. 


Table - 1 

Bayes point estimate of reliability function for m=n 
and varying mean failure rate 'a' and fixed T - K 
S.No. A R*sp(t) 


=n=3, 


failure rate 

'a' and fixed T - 

100 units. 

R'sp(t) 

R'ps(t) 

R'km(t) 

0.493 

O.lOO 

0.285 

0.341 

0.044 

0.178 

0.221 

0.019 

0.107 

0.136 

0.007 

0.063 

0.080 

0.003 

0.036 



Table - II 

Bayes point estimate of availability for m=n=5, r=k-j, s-2 
and varying mean failure rate ’a’ and fixed repair rate 'b' 


S.No. 

a 

b 

A* 

sp 

A* ^ 

ps 

A* 

^ km 

1. 

1 

1 

; 0.955 

0.1654 

0.392 

2. 

2 

1 1 

0.910 

0.105 

0.3-11 

3. 

3 

1 

0.888 

0.069 

0.249 

A 

4 

1 

i 

0.819 

0.046 

0.203 



5.4.4 CONCLUDING REMARKS 

In many practical situations, it is important to estimate the reliability 
and availability of the static system models when some operational 


experiences with the complete system are available as experimental data. 


: truncated failure and repair information have 


In the present chapter, time 
been used to update the analysis regarding the reliability and availability 


of such system. 

;-I reveal that the Bayes estimates ot R^Jt) and Rj,„/^) 


Table- 


;reases. Also tor fixed 


tend to be lower as the mean failure rate a mcr 
mean failure rate 'b' the Table II shows that the Bayes estimates of 

and A tend to be lower as mean 
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'P 


lean failure rate 'a' increases, rate a 


km 

increases. 


5.5.0PTIMIZATI0N OF SYSTEM AVAILABILITY 

IN THE BAYESIAN FRAMEWORK 

5.5.1. Introduction 

Optimization of reliability characteristics subjectto constraints on 
rbe number of components in static system models has been dre main 
concern of reliability engineers. Reliability and availability aie two imp 
characteristics of a system or device. For a system, availability 
as the ratio of mean time between failure (MTBF) to the mean time ' 
between failure (MTBF) plus mean time to repair (MTTR). Thei etoi e, 
both the reliability and maintainability aspects are considered for defining 
availability. Most of the reliability characteristics of various l.tetime 
distribution have been analysed in respect of their- statistical properties. 
However, inferences on the system availability are not many in literature. 

G..ayandLewis(l%7),MastersandLeivis[.^B7,, and Masters, Lewis 

and Kolarik ( 1992 ) have touched this aspect to some extent. Ihey 

devdlopedconfidenceinteiwalforavailabilityafterseparatelyestablishing 

advancement in technology in present time, it is reasonable to assume 
mndom variations in the parameters of failure time and repair time 

distributionduetoenvironme„talconditions.Thus,theBayesinaana,ysis 
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of the system availability become important which combines past 
experience with experimental data. In view of the above, the present 
section deals with the development of Bayesian tolerance limits for the 
evdilability of some basic system models, Subject to cosy constraints, 
these limits may be analysed to obtain tile optimum number ot components 
or subsystems in a system models with a preassigned confidence 

coefficientTimetruncatedsampleinformationhavebeenusedinanalysis. 


5 . 5 . 2 . Bayesian tolerance limits for system availability 

The Bayesian tolerance limits for system availability (A.) a~needed 

toattainatleastpavailability.withatleastconfidencecoefb^ 

be expressed as ; 

5.27) and the results in sections 5.2. the Bayesian 


Islow, using (5.27) 

admissible standards expressed in terms ot p 
osed to find the optimum number of components 

di fferent system models. 


j'oand y, these limits may be 


or subsystems 


in 
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5.5.3. (I) Series-parallel system 

For a series-parallel system, at least (100)/o Bayesia 
tolerance-limits for A, in (5.6) can be obtained from (5.27), i.e. using: 
(5.6) in (5.27), one gets, 

p[{i-(i-A)™}'’>p]^y 

Or P[A>{l-(l-p''" ^7 
Or P [ A <{ 1 -( 1 - p '''')'''”}]-^''7 

Now. on usingthe posterior distribution of A in (5.20) we have 




iAS+‘>-i(j_^)t«-'dA < 1-y 


.(5.28) 


B(s + b,r + a) 

Or l 7 __(s + b,r + a) <l-y 

.= ,-(, and l„(s + b,, «) is the well known 


Where z 

incomplete beta ratio. 

■ For attaining at 

(5.28) can be used to optimize m(n) with known n(m). 


least p availability with admissible y , the inequality 


5.5.4. (II) Parallel-series system 

Similarly on using (5.1 1 ) in (5.27), one gets, 
p[!l-(l-A-rihp]27 
Or p[Ah{i-ii-p) 


\ J 


; 

I, 


1 fa 
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. Or p[A<{i-(i-py''"r]^^-T 

, (i-(i-p)""'r 

B(s + b,r + a) •’ 

Or lz,(s + b,r + a)<l-y 
Inequality (5.29) can also be analysed as in the series-paiallel case. 


.(5.29) 


5.5.5. (HI) K-out-of-m System 

On substituting in, that limits for this system becomes 
arbitrarily chosen values of p and g have been summarized m table III 
and IV. These values can be analysed to optimize the values ot m and n 
for prc-assigned values of p andy in a series, parallel, series-parallel and 
parallel-series system models. Further, it is to be noted here that the 
,esults for the different system models in talbe 111 and IV show that an 
additional component in series adversely affects the values of p and 
However, these values attain a reasonable standard when an additional 



component is placed in a pai 


rallel system models. SimilaiW, the 


optimization of availability in the case 

adrieved by analysing the relevant inequalities. 


of a k-out-of-m system can be 
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Table III 

Bayesian tolerance limits for availability of a series and parallel system 


Series System 


Series System 


P 

Y 

n 

P 

Y 

N 

0.18 

0.34 

2 

0.60 

0.60 

o 

J 

0.45 

0.34 

1 

0.80 

0.60 

4 

0.10 

0.40 

2 

0.90 

0.60 

6 

0.10 

0.20 

3 

0.70 

r 0.50 

o 

J 

0.20 

0.10 

3 

0.70 

0.85 

5 

0.20 

0.20 

2 

0.70 

. 0.90 

6 


Table IV 

Bayesian tolerance limits for availability of a series - parallel and 

Parallel-series system 


Series-Parallel System 

Parallel-Series System 

|i 

P 

Y 

m 

n 

P 

y 

m 

n 

,1' 

1 

n i 


0.60 


9 

0.50 

0.50 

5 


0.40 

3 





; 

'f 

2 


0.60 

o 

o 

0.70 

0.50 

8 

0.30 

J 

9 




0.50 

0.40 

4 

5 

0.60 

0.54 

7 

2 

0.50 

0.50 

4 

4 

0.60 

0.60 

8 

o 

0.54 

0.70 

4 

2 

0.40 

0.27 

10 

4 

0.35 

0.70 

4 

3 

0.40 

0.40 

16 

4 
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